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feedbacks  in  the  region.  In  particular,  comparisons  of  RGPS  data  and  sea  ice  results  from 
ice-ocean  and  fully  coupled  regional  climate  models  have  been  made  to  evaluate  model 
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I.  INTRODUCTION 


A.  BACKGROUND  AND  MOTIVATION  FOR  SEA  ICE  RESEARCH 

The  Earth’s  climate  system  is  a  constantly  changing  environment,  which  can  be 
considered  in  balance  only  through  the  interdependence  of  its  dynamic  constituents. 
While  the  atmosphere  tends  to  capture  our  attention  on  a  daily  basis  through  changing 
weather  patterns,  it  is  the  ocean’s  variability  and  its  interaction  with  the  atmosphere  that 
provide  the  main  controls  on  large-scale  climate  patterns  (Bindoff  et  al.  2007). 
Understanding  air-sea  interactions  has  been  integral  to  our  comprehension  of  how  the 
climate  reacts  to  the  global  changes  that  have  been  observed  in  recent  history.  One  of  the 
emerging  areas  for  climate  research  is  to  study  the  Arctic,  where  a  layer  of  sea  ice  has 
perennially  altered  the  air-ocean  interface.  The  sea  ice  acts  as  a  barrier  between  the 
ocean  and  atmosphere,  and  the  kinematics  of  this  ice  greatly  impacts  the  heat  exchange, 
momentum  flux,  and  mass  balance  of  sea  ice  in  the  Arctic  (Rampal  et  al.  2009;  Moritz  et 
al.  2002).  Changes  to  sea  ice  coverage  due  to  kinematic  forcing  have  large  implications 
for  these  balances  and  climate  at  hemispheric  and  global  scales  more  generally  (Kwok 
and  Sulsky  2010;  Alexeev  2010). 

Over  the  last  several  decades,  observations  of  the  Arctic  point  towards  a  greatly 
diminished  ice  cap  due  to  a  variety  of  environmental  factors.  While  global  temperatures 
have  been  on  the  rise  over  the  last  century,  the  Arctic  atmosphere  has  warmed  by  almost 
two  times  the  global  average  in  the  last  100  years  (IPCC  2007a).  This  trend  has  been 
largely  attributed  to  changes  in  the  last  20  to  40  years,  where  surface  air  temperature 
increased  significantly  greater  than  previous  rates  (Overland  et  al.  2004;  Chylek  et  al. 
2009).  This  enhanced  temperature  increase  in  the  Arctic  has  been  well  documented  and 
has  been  attributed  to  a  positive  feedback  loop  unique  to  ice  covered  regions  (Perovich  et 
al.  2007;  Serreze  et  al.  2009).  As  the  climate  warms,  the  length  and  intensity  of  the  melt 
season  is  increased  which  leads  to  less  sea  ice  surviving  to  the  next  fall.  This  increases 
the  absorption  of  solar  energy  by  the  open  ocean,  which  increases  its  sensible  heat 
content  (Serreze  et  al.  2009).  This  additional  oceanic  heat  delays  the  formation  of  new 
ice  in  the  fall  and  winter,  promoting  enhanced  upward  heat  fluxes  into  the  atmosphere, 
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which  is  maintained  in  the  lower  troposphere  by  strong  low-level  stability  (Serreze  et  al. 

2009) .  This  polar  amplification  is  caused  in  part  by  the  dark  ocean  surface  left  behind 
from  the  melting  process,  which  has  a  much  lower  albedo  (-0.07),  a  measure  of  surface 
reflectivity  as  a  fraction  of  the  incident  solar  radiation,  than  bare  sea  ice  or  snow-covered 
ice  (0.65  to  0.85)  (Perovich  et  al.  2002).  Larger  open  water  surface  and  heat  content  of 
the  upper  ocean  can  not  only  further  melt  sea  ice,  but  also  alter  oceanic  and  atmospheric 
circulation  patterns,  and  add  more  water  vapor  into  the  lower  atmosphere  by  releasing 
heat  stored  in  the  mixed  layer  before  freezing  in  the  fall.  This  amplifies  the  impact  on  air 
temperature  increase  as  water  vapor  causes  a  regional  greenhouse  effect  to  the  lower 
atmosphere  (Langen  and  Alexeev  2007).  The  positive  feedback  loop  has  helped  to 
perpetuate  the  increased  warming  conditions  recently  observed  in  the  Arctic,  and  is  one 
of  the  reasons  scientists  fear  that  this  region  might  be  approaching  a  new  regime  of 
seasonal  sea  ice  coverage. 

The  primary  observational  evidence  of  this  negative  trend  in  sea  ice  has  been  in 
the  analysis  of  sea  ice  extent  from  satellite  passive  microwave  imagery  since  the  late 
1970s  (Johannessen  et  al.  2004;  Comiso  et  al.  2008;  Kwok  et  al.  2008;  Kwok  and  Sulsky 

2010) .  Sea  ice  extent  is  defined  as  the  area  of  ocean  with  at  least  15  percent  sea  ice 
coverage  (within  a  satellite  footprint).  This  definition  allows  for  the  inclusion  of  melt 
ponds  on  the  ice  pack  to  alleviate  the  inability  of  satellite  imagery  to  differentiate 
between  water  from  these  features  (with  sea  ice  below)  and  that  of  open  ocean.  Based  on 
this  definition,  sea  ice  in  the  Arctic  region  has  extended  from  14  to  16  million  square 
kilometers  in  late  winter  (March)  and  from  six  to  eight  million  square  kilometers  by  the 
end  of  each  summer  (September)  during  1979-2000  (National  Snow  and  Ice  Data  Center 
(NSIDC)  2011).  In  the  last  decade  however,  sea  ice  extent  in  September  has  been 
trending  between  4  to  6  million  square  kilometers  of  coverage  (Figure  1).  This  decline  in 
ice  extent  was  one  of  the  first  indicators  of  accelerating  climate  warming  in  the  Arctic, 
which  can  be  linked  to  the  positive  feedback  loop  related  to  the  ice-albedo  effect. 
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Figure  1.  Arctic  sea  ice  extent  through  September  2011  (From  National  Snow 
and  Ice  Data  Center  [NSIDC]  2011) 

1.  Atmospheric  Influences 

Sea  ice  loss  is  not  limited  to  reduction  in  just  areal  extent,  however,  as  export  of 
old  and  thick  ice  out  of  the  Arctic  basin  and  thennodynamic  forcing  from  the  ocean  and 
atmosphere  has  led  to  a  decrease  in  overall  ice  thickness  as  well.  The  thickness  of  sea  ice 
influences  characteristics  of  the  ice  cover  as  a  thermodynamic  barrier  as  well  as  greatly 
impacting  the  dynamic  nature  of  ice  motion.  Dynamic  processes  in  sea  ice  are 
determined  by  ice  thickness,  strength,  and  the  coupling  of  forces  from  the  air  above  and 
ocean  below.  Ice  motion  has  been  linked  to  large-scale  atmospheric  circulation  patterns 
associated  with  the  Arctic  Oscillation  (AO)  or  the  North  Atlantic  Oscillation  (NAO),  also 
referred  to  as  the  Northern  Annular  Mode  (NAM)  (Thompson  and  Wallace  1998;  Serreze 
et  al.  2007).  The  AO  develops  as  atmospheric  pressure  patterns  around  Iceland  and 
Azores  dominate  the  regional  flow  field  for  extended  periods  of  time.  A  positive  phase 
of  the  AO  occurs  when  a  negative  sea  level  pressure  (SLP)  anomaly  remains  over  the 
pole  with  a  positive  SLP  anomaly  dominating  the  mid-latitudes,  and  a  negative  phase 
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occurs  when  this  pattern  is  weakened  or  reversed.  The  NAO  is  seen  as  a  pattern  on  the 
Atlantic  side  that  works  to  amplify  or  diminish  the  effects  of  the  AO  through  the 
covariability  between  the  strength  of  the  Icelandic  Low  and  the  Azores  High  (Serreze  et 
al.  2007).  During  a  positive  NAO,  a  strong  low  pressure  system  is  present  in  the  North 
Atlantic  near  Iceland,  and  a  strong  high  pressure  system  is  seen  in  the  subtropics  near  the 
Azores,  which  combined  increase  the  intensity  of  the  Westerlies  across  the  midlatitudes 
deflecting  the  stonn  track  further  north  into  the  Arctic  (Maslanik  et  al.  1996;  McPhee  et 
al.  1998).  Changes  in  storm  tracks  can  either  increase  (+NAO)  or  decrease  (-NAO)  the 
energy  and  atmospheric  mass  imported  to  the  Arctic,  which  determines  the  overall 
strength  of  the  AO  (Serreze  et  al.  2007). 

While  the  AO/NAO/NAM  have  an  impact  on  weather  patterns  and  storm 
propagation  in  the  mid-latitudes,  in  the  Arctic  they  also  impact  the  flow  of  sea  ice  as 
pressure  differences  can  change  prevailing  wind  patterns.  During  a  negative  phase  of  the 
AO,  a  positive  pressure  sets  up  an  anticyclonic  (clockwise),  zonal  flow  in  the  Beaufort 
Gyre,  which  is  believed  to  cause  the  ice  to  converge  and  raft  in  this  region  due  to  Coriolis 
forcing.  This  process  would  allow  sea  ice  to  increase  in  thickness  and  remain  longer  in 
the  Arctic  basin,  which  would  lead  to  an  increase  in  the  presence  of  older,  thicker  ice. 
During  a  positive  phase  of  the  AO,  this  atmospheric  pattern  is  weakened  or  reversed 
which  would  tend  to  allow  ice  to  flow  more  freely  along  the  principle  pathway  out  of  the 
Arctic  through  Fram  Strait,  between  Greenland  and  Svalbard  (Figure  2).  From  the  late 
1980s  through  the  mid-1990s,  the  AO  entered  a  strongly  positive  phase,  which  was  seen 
as  the  primary  atmospheric  forcing  that  lead  to  an  increased  ice  export  through  Fram 
Strait  and  the  decrease  in  sea  ice  thickness  as  it  was  thought  to  induce  both  the  warming 
and  export  necessary  to  account  for  the  observed  trends  (Thompson  and  Wallace  1998; 
Rigor  and  Wallace  2004).  Since  2002,  however,  the  AO  has  shifted  back  to  a  largely 
neutral  to  negative  phase,  while  sea  ice  extent  record  lows  have  continued  to  be  observed 
throughout  the  last  decade  (Maslanik  et  al.  2007).  This  has  led  climate  scientists  to 
search  for  a  new  explanation  for  the  rapid  ice  loss  in  the  Arctic,  or  to  question  whether  or 
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not  sea  ice  presence  in  the  Arctic  has  reached  a  new  regime  of  the  overall  forcing  that 
determines  ice  motion  and  seasonal  survivability  of  ice  within  the  region  (Lindsay  and 
Zhang  2005). 


Figure  2.  Sketch  of  a  (a)  positive  phase  and  a  (b)  negative  phase  of  the  Arctic 
Oscillation  induced  atmospheric  patterns  (From  Wallace  2011). 

A  refined  look  at  the  atmospheric  forcing  has  revealed  that  the  second  leading 

mode  of  the  sea  level  pressure  (SLP)  anomalies  in  the  Arctic  winter  (Oct  -  Mar)  might  be 

gaining  importance  recently  with  regard  to  sea  ice  motion,  when  compared  to  the  AO 

(Wu  et  al.  2006).  This  pattern  is  known  as  the  Arctic  Dipole  Anomaly  (DA),  or  simply 

Arctic  Dipole  (AD)  in  other  publications,  which  is  determined  by  a  more  local  SLP 

fluctuation  where  pressure  centers  over  the  Laptev  Sea  (near  Siberia)  and  the  Greenland 

Sea  cycle  between  positive  and  negative  phases  in  a  dipole  structured  pattern  that  induces 

meridional  flow  (Watannabe  et  al.  2006;  Overland  and  Wang  2010).  While  the  DA  only 

accounts  for  14-19%  of  the  Arctic’s  variance  in  SLP  (as  opposed  to  59-63%  from  AO), 

it  remains  important  to  the  advection  of  sea  ice  within  and  out  of  the  Arctic  basin  since  it 

drives  the  position  of  these  local  pressure  centers.  A  positive  phase  of  the  DA  is 

associated  with  a  negative  SLP  anomaly  on  the  North  American  side  of  the  Arctic 

(Overland  and  Wang  2010).  While  the  anomalous  centers  of  this  pattern  vary  in  space, 

when  the  positive  SLP  anomaly  is  over  the  western  Arctic  and  the  negative  SLP  anomaly 
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is  over  the  Greenland  Sea,  atmospheric  forcing  increases  ice  flow  along  the  Transpolar 
Drift  Stream  (TDS)  due  to  the  meridional  nature  of  the  forcing  (Wu  et  al.  2006;  Richter- 
Menge  and  Overland  2009).  This  forcing  also  sets  up  an  acceleration  of  the  melting 
process  in  the  western  Arctic  as  the  oceanic  heat  flux  through  Bering  Strait  is  increased, 
causing  delayed  formation  of  new  ice  in  fall  and  earlier  onset  of  melting  in  the  Chukchi 
and  East  Siberian  Seas  in  spring  (Wang  et  al.  2009). 


62% EOF  1 
f 


{  <-  -I 


Figure  3.  Regression  maps  of  the  first  two  leading  modes  to  the  (a  and  b)  winter 
and  (c  and  d)  summer  mean  Northern  Hemisphere  SLP  field  using  the 
NCEP  Reanalysis  dataset  from  1948  to  2008.  Contour  intervals  are  0.5 
hPa  (see  color  bars).  The  black  arrows  in  (a)  and  (c)  indicate  the  cyclonic 
(anticyclonic)  wind  anomaly  during  the  +AO  (-AO).  In  (b)  and  (d),  the 
black  arrows  indicate  wind  anomalies  between  the  western  and  eastern 
Arctic  during  +DA  and  -DA  phases  that  accelerate  or  decelerate  the  TDS 
(in  red-dashed  arrows),  respectively 
(From  Wang  et  al.  2009). 

The  DA’s  relative  importance  to  sea  ice  export  and  melt  was  expanded  upon  by 
Wang  et  al.  in  2009,  as  they  ran  global  climate  models  (GCMs)  with  forcing  from  the 
National  Center  for  Environmental  Prediction  and  the  National  Center  for  Atmospheric 
Research  (NCEP/NCAR)  reanalysis  datasets  of  years  (1948-2004)  that  included  all 
potential  combinations  of  both  AO  and  DA.  They  found  that  regardless  of  the  sign  of  the 

AO,  sea  ice  export  was  greatly  increased  during  years  with  a  positive  phase  of  the  DA. 
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This  is  likely  due  to  the  directional  dependence  of  the  TDS,  which  parallels  the 
Greenwich  meridian  during  years  of  positive  DA,  leading  to  more  sea  ice  export  out  of 
the  Arctic  basin  due  to  the  strong  low  pressure  anomaly  in  the  Greenland  Sea.  During 
years  of  negative  DA,  more  ice  tends  to  be  held  within  the  basin  as  the  Beaufort  gyre  is 
strengthened  and  flow  along  the  TDS  is  weakened,  which  accumulates  more  ice  and  fresh 
water  in  that  region  (Wu  et  al.  2006).  A  negative  DA  pattern  can  however  push  more  ice 
back  into  the  Arctic  leaving  more  open  water  in  the  North  Atlantic,  which  can  increases 
the  amount  of  heat  storage  in  the  ocean  and  eventually  impacts  the  atmospheric 
circulation  through  release  of  this  excess  heat  in  fall,  preconditioning  ice  for  more  retreat 
in  following  seasons  (Overland  and  Wang  2010). 

2.  Sea  Ice  Thickness 

Increased  sea  ice  flow  out  of  the  Arctic  correlated  to  these  atmospheric  forcing 
patterns  in  recent  years  not  only  increases  the  ice-albedo  effect  through  an  increase  in  the 
area  of  open  ocean,  but  it  greatly  impacts  the  distribution  of  older  and  thicker  ice. 
Observations  of  average  annual  sea  ice  area  flux  from  satellite  imagery  from  1978-2002 
show  a  mean  rate  of  approximately  866,000  km  /year  southward  through  Fram  Strait 
(Kwok  et  al.  2004a).  This  average  rate  was  heavily  influenced  by  recent  increase  in  the 
fluxes  where  the  net  exported  areal  extent  was  -400,000  km  greater  in  the  1990s  than  in 
the  1980s  (Kwok  et  al.  2004a).  The  increased  export  of  almost  half  the  annual  areal 
average  has  influenced  the  thickness  and  age  distribution,  as  some  of  the  fonnerly 
resilient  ice  has  been  flushed  out  of  the  Arctic  during  this  period.  Data  derived  from 
upward  looking  sonar  (ULS)  thickness  distributions  across  Fram  Strait  over  an  eight  year 
period  from  1991  -  1998  provides  an  estimated  average  annual  ice  volume  flux  of  -2200 
km  of  sea  ice  through  Fram  Strait  (Kwok  et  al.  2004a).  This  volume  flux  corresponded 
to  a  period  of  overall  thickness  decrease  of  0.45m  in  the  Arctic  region  through  the  early 
2000s,  with  a  greater  decrease  of  0.55m  in  multi-year  (MY)  ice  when  compared  to 
previous  annual  averages  (Kwok  et  al.  2004a).  This  increased  export  and  decrease  in  sea 
ice  thickness  occurred  over  a  period  of  consistently  high  positive  AO  in  both  fall  (Oct- 
Dee,  1988-1994)  and  winter  (Jan-Mar,  1989-1997),  which  helped  to  trigger  the  observed 

sea  ice  response  (Thompson  and  Wallace  1999;  Overland  and  Wang  2010). 
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Ice  thickness  is  difficult  to  observe  at  basin  scale  as  transects  of  submarine 
observations  from  under  the  ice  cap  are  limited  snapshots  in  time,  moored  buoys  are  only 
possible  at  particular  points  in  space,  and  ice  cores  are  limited  in  both  space  and  time. 
The  Scientific  Ice  Expeditions  (SCICEX)  program,  started  in  the  1990s  by  the  U.S.  Navy, 
sought  to  address  the  need  for  observational  data  by  involving  scientists  into  the  planning 
and  actual  measurement  of  ice  thickness  during  dedicated  science  cruises.  While  still 
limited  in  basin  wide  coverage,  this  data  collection  program  provided  the  best  ice 
thickness  data  set  to  date.  It  has  allowed  for  comparison  with  previous  data  released  by 
the  Navy  from  cruises  that  were  declassified  from  the  Cold  War  era  operations  from  the 
1950s  to  the  late  1970s  (SCICEX  Science  Advisory  Committee  2010).  Analysis  of  this 
data  revealed  that  the  ice  draft  along  similar  transects  from  the  central  Arctic  basin  had 
decreased  by  over  a  meter  in  the  1990s  from  the  over  three  meter  averages  taken  from 
submarine  cruises  in  previous  decades  (Rothrock  et  al.  1999).  This  significant  reduction 
in  thickness  corresponded  to  a  total  sea  ice  volume  reduction  of  around  40%  in  a  matter 
of  a  few  decades,  which  was  largely  attributed  to  large  heat  fluxes  from  the  ocean  and 
atmosphere  (Rothrock  et  al.  2003).  An  alternative  explanation  of  this  volume  reduction 
(Holloway  and  Sou  2002)  has  been  that  the  measured  thickness  changes  were  due  to  the 
redistribution  of  ice  outside  the  submarine  tracks  toward  the  Canadian  Archipelago  and 
Siberian  shelves.  This  shows  that  while  the  data  collected  during  the  1990s  are  a  vast 
improvement  over  the  sporadic  data  available  from  previous  naval  operations,  they  are 
still  limited  spatially  by  narrow  range  transects  of  observations  and  temporally  due  to  the 
exercise  sampling  only  in  late  summer  months. 

To  combat  these  spatial  and  temporal  limitations  in  ice  thickness  data,  techniques 
have  been  developed  recently  to  estimate  sea  ice  age  and  thickness  through  a 
measurement  of  freeboard  from  satellite  imagery  taken  from  the  Ice,  Cloud,  and  land 
Elevation  Satellite  (ICESat)  campaign,  which  can  cover  nearly  the  entire  basin  (up  to 
84°N)  and  can  be  measured  throughout  the  year  (Kwok  et  al.  2004b;  2007).  Freeboard  is 
the  height  of  ice  and  snow  that  is  above  mean  sea  level  (MSL),  from  which  the  ice 
thickness  can  be  estimated.  This  approach  assumes  hydrostatic  equilibrium  and  that  snow 
depth  and  ice  density  are  known  (Kwok  and  Cunningham  2008).  Snow  characteristics 
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are  difficult  to  differentiate  from  satellite  imagery  alone,  so  initial  estimates  have  to  be 
taken  from  climatology  and  meteorological  products  to  estimate  snow  density  and 
distribution  throughout  a  given  season.  While  limitations  in  estimating  the  seasonal  snow 
cycle  remain,  ICESat  analysis  using  these  assumptions  has  given  rise  to  regional 
estimates  with  an  uncertainty  in  thickness  calculations  of  ~0.7m.  Recent  ICESat 
observations  and  subsequent  estimates  of  ice  thickness  have  been  validated  by 
comparison  to  submarine  and  buoy  data.  The  variance  between  ICESat  draft  and 
submarine  draft  was  calculated  to  be  0.42m  (Kwok  and  Rothrock  2009),  while 
comparisons  with  buoy  data  estimates  are  within  0.5m  of  buoy  observations  (Kwok  and 
Cunningham  2008).  When  used  in  conjunction  with  these  datasets,  satellite  derived  ice 
thickness  offer  an  overall  picture  of  basin-wide  characteristics  of  the  Arctic  sea  ice 
thickness  distribution  and  trends  (Figure  4)  (Kwok  et  al.  2009;  Kwok  and  Rothrock 
2009).  Additional  satellite  imagery  techniques  will  be  further  described  in  following 
chapters. 


1970  1980  1990  2000  2010 


Figure  4.  (a)  Winter  Arctic  Ocean  sea  ice  thickness  from  ICESat  (2004-2008). 

(b)  Interannual  changes  in  ice  thickness  from  the  submarine  (RA)  and 
ICESat  campaigns  (After  Kwok  and  Rothrock  2009;  Kwok  and  Sulsky 
2010). 
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3.  Sea  Ice  Age 

Given  the  difficulty  in  measuring  sea  ice  thickness,  ice  age  is  used  as  a  proxy  for 
thickness  estimations.  As  blocks  of  ice  remain  in  seawater  over  several  seasons,  pockets 
of  brine  created  during  fonnation  get  flushed  out  reducing  the  salt  content  and  increasing 
the  strength  of  the  remaining  sea  ice  (Tucker  et  al.  1992).  Due  to  a  different  salt  content, 
differences  in  the  brightness  temperature  (a  measure  of  the  emissivity  of  the  surface 
material)  of  first  year  (FY)  and  multi-year  (MY)  ice  (ice  that  lasts  through  at  least  one 
melt  season)  can  be  detected  and  tracked  by  satellite  images  from  Scanning  Multichannel 
Microwave  Radiometer  (SMMR),  Special  Sensor  Microwave/Imager  (SSM/I),  and  the 
series  of  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  sensors  (Maslanik  et  al. 
2007).  Ice  age  can  be  derived  from  a  combination  of  satellite  imagery  and  validated  with 
ocean  buoy  data  where  Lagrangian  grid  cells  are  tracked  throughout  the  Arctic  Ocean  and 
considered  to  maintain  ice  through  a  melt  season  as  long  as  they  contain  at  least  15%  sea 
ice  concentration  (Fowler  et  al.  2004).  These  cells  when  tracked  over  a  number  of  years 
detennine  ice  age  within  a  given  region.  This  parameter  is  important  to  understanding 
sea  ice  dynamics  as  ice  strength  generally  increases  with  thickness  and  age  due  to  brine 
rejection  from  within  the  ice  (Kovacs  1996).  This  process  not  only  impacts  the  salinity  of 
the  ocean  and  the  fresh  water  storage  within  the  Arctic,  but  also  greatly  impacts  the 
rigidity  and  mobility  of  the  pack. 

Observations  have  shown  that  since  the  1980s,  MY  ice  has  been  diminished  from 
covering  -80%  of  the  Arctic  in  that  decade  to  only  -30%  of  the  area  since  the  1990s 
(Rigor  and  Wallace  2004,  Drobot  et  al.  2008,  Kwok  and  Cunningham  2010).  Along  with 
sea  ice  export  and  melt,  redistribution  of  thick,  MY  ice  has  shifted  the  ice-covered  area  in 
the  Arctic  toward  younger  and  thinner  ice,  which  is  more  prone  to  move  around,  ridge, 
and  melt  during  the  next  summer.  This  younger  and  thinner  ice  is  also  more  likely  to 
experience  abrupt  changes  in  ice  export  under  the  same  atmospheric  and  oceanic  forcing 
that  it  could  resist  in  the  past  (Maslanik  et  al.  2007).  Since  the  all-time  record  minimum 
in  sea  ice  extent  in  2007,  areas  such  as  the  central  Arctic  and  the  Canadian  Archipelago 
that  were  relatively  stable  in  the  past  have  now  lost  vast  amounts  of  old  and  thick  ice, 
which  points  to  an  overall  shift  in  the  Arctic  sea  ice  cover  (Figure  5)  (Maslanik  et  al. 
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2011).  A  significant  question  remaining  is  whether  sea  ice  dynamics — specifically,  the 
convergence  of  thinner  MY  ice  coverage — could  be  misconstrued  as  a  net  loss  of  ice 
from  melt  or  export,  while  it  actually  represents  a  reduction  in  the  ice  area  (Kwok  and 
Cunningham  2010). 


Figure  5.  Median  age  of  March  sea  ice  in  1985  (left)  and  2011  (right).  Overall,  the 
proportion  of  young  ice  (light  blue)  has  increased.  By  March  2011,  ice 
more  than  four  years  old  (dark  blue)  accounts  for  less  than  10  percent  of 
the  Arctic  ice  cover  (From  Maslanik  et  al.  2007) 


B.  INTENT  OF  STUDY 

Sea  ice  extent,  age,  and  thickness  have  been  on  the  decline  in  the  Arctic  over  the 
past  several  decades,  and  record  low  ice  concentration  and  thickness  are  becoming 
commonplace  at  the  end  of  each  summer.  Data  from  satellites,  buoys,  ice-tethered 
profilers  (ITPs),  and  submarine-based  sensors  all  point  to  an  increasing  trend  toward 
seasonal  ice  coverage  in  the  region.  Understanding  the  mechanisms  that  drive  these 
trends  is  the  key  to  determining  whether  this  pattern  will  continue,  or  if  natural  variability 
is  capable  of  forcing  a  reversal.  There  are  two  types  of  forces  that  impact  sea  ice 
response  to  climatic  changes  in  the  Arctic:  thermodynamic  and  dynamic  forcing. 
Thermodynamic  processes  involve  changes  in  surface  air  temperature  (SAT),  radiative 
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fluxes,  and  oceanic  heat  transport,  while  dynamic  processes  involve  changes  in  sea  ice 
motion  and  circulation  in  response  to  oceanic  and  atmospheric  forcing  (Serreze  et  al. 
2007).  This  study  focuses  on  the  kinematics  of  sea  ice  in  the  Arctic  to  help  determine 
model  capabilities  in  representing  the  behavior  of  sea  ice. 

While  few  contemporary  studies  have  analyzed  these  processes,  the  scientific 
community  has  recently  begun  to  look  at  the  dynamical  forcing  that  detennines  the 
distribution  and  mechanical  characteristics  of  sea  ice  (Rothrock  and  Zhang  2005;  Kwok 
et  al.  2008;  Rampal  et  al.  2009).  Understanding  the  relative  contribution  of  dynamic 
processes  to  sea  ice  state  will  help  improve  the  current  understanding  of  the  role  sea  ice 
plays  in  the  Arctic  response  to  climate  change.  Advances  in  computing  power  have 
offered  the  possibility  of  modeling  of  fine  scale  processes  and  forcing  to  better  represent 
the  sea  ice  thickness  distribution  in  the  Arctic  and  the  survivability  of  sea  ice  in  a  given 
year.  Through  the  analyses  of  observed  and  modeled  dynamic  processes  of  sea  ice  drift 
and  deformation  presented  here,  we  aim  to  improve  the  understanding  of  the  recent 
changes  observed  in  Arctic  sea  ice. 

Numerical  modeling  remains  a  key  method  of  analysis  and  synthesis  of  air-ocean- 
sea  ice  interactions.  While  the  preponderance  of  observations  remain  to  be  of 
atmospheric  processes  and  model  re-analyses,  which  have  led  to  many  studies  on  the 
atmospheric  forcing  of  Arctic  sea  ice  (Serreze  et  al.  2003,  2007;  Stroeve  et  al.  2005, 
2011;  Perovich  et  al.  2007),  a  better  understanding  of  the  sea  ice  dynamical  response  to 
this  forcing,  and  its  influence  on  thickness  distribution  has  been  largely  overlooked. 
When  looking  at  the  overall  heat  budget  and  the  increased  energy  required  to  melt  or 
deform  Arctic  sea  ice  in  the  manner  that  has  been  recently  observed,  it  has  become  clear 
that  atmospheric  forcing  alone  is  not  sufficient  to  account  for  this  loss  (Francis  et  al. 
2005;  Maslowski  et  al.  2007).  Sea  ice  deformation  and  oceanic  forcing  therefore  have  to 
be  more  closely  analyzed  to  determine  their  role  in  the  recent  decline  of  sea  ice  extent 
and  thickness. 
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C.  OVERVIEW 

This  study  analyzes  model  results  from  several  regional  Arctic  models  developed 
and  maintained  at  the  Naval  Postgraduate  School  (NPS).  The  model  domain  includes  the 
sub-Arctic  North  Pacific  and  North  Atlantic  Oceans,  the  Arctic  Ocean,  the  Canadian 
Archipelago  and  the  Nordic  Seas.  The  main  focus  is  on  the  effects  of  dynamic  forcing  on 
sea  ice  thickness,  distribution,  and  variability  within  the  Arctic  region.  Seasonal, 
regional,  and  interannual  variability  in  sea  ice  motion  fields  and  deformation  rates  are 
investigated.  The  model  results  for  the  Arctic  region  are  synthesized  with  observational 
data  from  RADARSAT  imagery  to  make  comparisons  to  observed  sea  ice  deformation 
based  on  previous  studies. 

This  thesis  is  organized  into  the  following  chapters:  Chapter  II  describes  U.S. 
Navy  relevance  of  sea  ice  research  in  the  Arctic;  Chapter  III  describes  the  models  and 
data  analysis  methodology  used  to  conduct  this  research;  Chapter  IV  presents  results  and 
provides  discussion;  Chapter  V  consists  of  a  summary  and  conclusions  and;  Chapter  VI 
makes  recommendations  for  future  research. 
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II.  NAVY  RELEVANCE 


International  interest  in  the  state  of  the  Arctic  has  grown  in  recent  years  as 
observations  of  increasingly  open  ocean  area  in  summertime  have  validated  what  had 
been  reported  for  decades  by  the  Intergovernmental  Panel  on  Climate  Change  (IPCC). 
The  IPCC  consists  of  an  international  scientific  body  established  by  the  United  Nations 
Environment  Programme  (UNEP)  and  the  World  Meteorological  Organization  (WMO)  in 
1988  to  report  on  the  current  state  of  knowledge  on  climate  change  and  its  potential 
environmental  and  socioeconomic  impacts  (IPCC  2012).  This  group  is  an  open 
consortium  of  scientists  that  compiles  a  report  periodically  with  the  intent  of  building 
consensus  on  drivers  of  observed  climate  change  in  a  peer-reviewed  environment.  With 
the  release  of  its  fourth  assessment  report  in  2007  (AR4),  the  IPCC  gained  widespread 
geopolitical  attention  with  its  assertion  that  anthropogenic  forcing,  specifically  the 
extensive  release  of  greenhouse  gases  (GHG)  from  human  activities,  was  causing  an 
amplified  impact  on  melting  ice  and  warming  temperatures  in  the  Arctic  region  (IPCC 
2007a). 

Although  the  United  States  (U.S.)  has  territorial  claims  to  resources  extending 
from  its  coast  and  has  operated  in  the  Arctic  for  exploratory,  scientific,  and  military 
purposes,  it  has  lacked  clearly  defined  guidance  for  the  management  of  its  interests  in  the 
region  until  recently.  Since  the  Cold  War  concluded  in  the  early  1990s,  the  U.S.  interest 
in  the  Arctic  from  a  military  perspective  has  been  on  the  decline  as  the  former  Soviet 
Union  collapsed  and  the  emerging  Russian  Federation  has  developed  improved  relations 
with  the  U.S.  government.  Apart  from  the  U.S.  and  Russian  dynamic,  including  their 
geopolitical  history  and  close  proximity  in  the  Bering  Strait,  the  Arctic  maintains  a 
unique  set  of  geographic  challenges  based  on  the  number  of  nations  that  claim  national 
interest  in  the  region.  The  Arctic  Ocean  is  confined  by  remote  landmasses  of  five 
stakeholder  nations  including  the  U.S.,  Canada,  Russia,  Norway,  and  Denmark 
(Greenland)  (Figure  6).  These  countries  along  with  Finland,  Iceland,  Sweden,  and 
indigenous  representatives  developed  what  is  known  as  the  Arctic  Council  in  1996  as  a 
forum  for  international  cooperation  with  the  intent  of  protecting  national  interests  and 
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regulating  disputes  with  a  focus  on  sustainable  development  and  environmental 
protection  (Arctic  Council  1996).  This  body  is  limited  in  its  authority,  however,  to 
oversee  legitimate  international  agreements  or  impose  legally  binding  restrictions  on 
nations  for  abusing  natural  resources  or  disregarding  sovereignty  claims  due  to  the  U.S. 
refusal  to  accept  security  concerns  as  part  of  the  Council’s  charter  (Borgerson  2008). 


Figure  6.  Map  of  the  Arctic  Region  showing  the  Arctic  Ocean,  adjacent  seas  and 
political  boundaries.  Red  line  shows  average  10°C  isotherm  for  July 
(From  U.S.  Central  Intelligence  Agency,  2012). 
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A.  STRATEGIC  GUIDANCE 

The  leading  international  treaty  concerning  maritime  affairs  including  the  rights 
and  responsibilities  of  nations  for  their  use  of  the  ocean,  as  well  as  guidelines  for 
economic  and  environmental  resources  within  the  world’s  oceans  is  the  United  Nations 
Convention  on  the  Law  of  the  Sea  (UNCLOS).  This  treaty,  which  was  internationally 
recognized  in  1994,  defines  the  rights  of  navigation  through  national  and  international 
waters  and  sets  up  key  provisions  for  long-disputed  territorial  claims  including 
continental  shelf  extent  and  exclusive  economic  zones.  Over  the  last  century,  many 
nations  have  debated  the  definitions  that  detennine  their  rights  to  surrounding  waters,  and 
the  geographical  limits  of  these  regions  have  been  extended  multiple  times.  That  which 
was  once  a  three  nautical  mile  extent  from  the  coastal  reaches  of  a  nation  has  been 
pushed  out  to  a  200  nautical  mile  exclusive  economic  zone  in  which  mineral  resources 
and  fisheries  in  this  oceanic  region  can  be  exploited  and  protected  by  that  state.  While 
the  United  States  agrees  with  much  of  this  document,  Congress  has  failed  to  ratify 
UNCLOS  based  on  a  dispute  in  Part  XI  of  the  Convention,  which  relates  to  mineral 
definitions  beyond  the  twelve  mile  territorial  waters  of  a  state,  that  has  been  deemed  to  be 
in  conflict  with  American  economic  and  security  interests.  Based  on  this  discrepancy, 
the  U.S.  is  not  a  signatory  nation  to  this  unique  document  and  therefore  is  not  currently 
part  of  the  UN  commission  which  would  have  jurisdiction  over  debated  rights  to 
resources  in  regions  such  as  the  Arctic  (Borgerson  2008). 

Despite  this  discord,  recent  guidance  has  been  initiated  due  to  gaining  acceptance 
of  the  IPCC  report,  growing  international  recognition  for  the  validity  of  its  contents,  and 
the  increased  interest  in  the  utilization  of  the  vast  resources  in  the  Arctic.  In  2009, 
President  George  W.  Bush  signed  the  National  Security  Presidential  Directive  (NSPD)  66 
/  Homeland  Security  Presidential  Directive  (HSPD)  25.  This  directive  was  the  first  in  a 
series  of  strategic  guidance  documents  that  began  to  define  the  national  security  interests 
and  objectives  for  U.S.  diplomacy  and  potential  operational  capabilities  in  the  Arctic. 
Along  with  its  support  to  ratify  UNCLOS,  this  mandate  established  requirements  for  the 
U.S.  to  develop  and  maintain  greater  capacity  to  protect  sovereign  interests  extending 
from  our  boarders  to  the  extent  of  the  continental  shelf,  and  to  increase  maritime  domain 
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awareness  and  global  mobility  through  an  enhanced  ability  to  operate  within  the  Arctic 
(Task  Force  Climate  Change  /  Oceanographer  of  the  Navy  2011). 

In  an  effort  to  address  the  impact  of  the  environmental  concerns  outlined  in  the 
Presidential  Directives  on  naval  operations,  the  Chief  of  Naval  Operations  (CNO), 
Admiral  Gary  Roughead,  created  Task  Force  Climate  Change  (TFCC)  in  May  2009. 
This  task  force  was  placed  under  the  direction  of  the  Oceanographer  of  the  Navy,  who 
was  charged  with  an  initial  goal  to  develop  and  implement  strategic  policy  initiatives  for 
the  Arctic  specifically,  and  then  to  follow  up  with  a  roadmap  for  Department  of  Defense 
global  climate  change  responses  more  generally  (TFCC  /  Oceanographer  of  the  Navy 
2009).  The  first  actionable  document  produced  was  the  Navy  Arctic  Roadmap,  which 
delineated  a  list  of  action  items  to  address  both  the  strategic  and  policy  implications  of 
changes  in  the  Arctic,  as  well  as  provided  a  vehicle  for  the  Navy  to  assess  and  predict 
future  environmental  changes  in  the  region.  From  a  strategic  perspective,  this  document 
focused  on  determining  the  potential  impact  of  increased  shipping  activity  from  renewed 
interest  in  natural  resource  exploration,  research  initiatives,  fishing  migrations  north,  and 
eco-tourism  into  a  more  environmentally  approachable  Arctic.  Increased  traffic  from 
these  activities  into  this  still  harsh  environment  along  with  the  recent  opening  of 
interconnected  sea  routes  of  the  Northwest  Passage  and  the  Northern  Sea  Route  across 
the  Arctic  in  2008  has  rightfully  alarmed  those  who  have  been  operating  in  the  region  for 
some  time  (Figure  7).  Increased  activity  will  undoubtedly  bring  with  it  the  potential  for 
environmental  disasters  such  as  oil  spills  and  increased  search  and  rescue  operations  as 
unprepared  ships  encounter  the  fast-moving  ice  in  the  marginal  ice  zone  (MIZ).  This, 
along  with  the  lack  of  accurate  mapping  in  the  region,  due  to  extensive  ice  coverage  in 
the  past,  will  likely  lead  to  the  need  for  improved  cooperation  between  Arctic  States  in 
order  to  effectively  support  the  increased  operational  safety  response  requirements. 
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Figure  7.  Potential  Summer  Shipping  Lanes  (From  Arctic  Marine  Shipping 
Assessment  (AMSA),  2009) 

B.  ARCTIC  NATIONS 

Stakeholders  in  the  Arctic  all  have  complex  individual  and  collective  relationships  with 

their  Arctic  neighbors.  The  Navy  Arctic  Roadmap  listed  the  development  of  strong 

cooperative  interagency  and  international  partnerships  as  one  of  the  primary  goals  of  this 

renewed  focus  on  the  Arctic.  The  importance  of  strengthening  these  relationships  is 

motivated  by  recent  extension  claims  and  disputes  over  territorial  rights  of  the  Arctic 

region.  The  Russians  have  been  attempting  for  over  a  decade  to  claim  that  the 

underwater  Lomonosov  Ridge,  which  extends  across  the  Arctic  basin,  is  an  extension  of 

their  continental  margin  and  the  resources  it  contains  should  therefore  be  considered 

exclusive  property  (O’Rourke  2011).  The  annex  of  this  undersea  area  would  grant  Russia 

over  50%  of  the  Arctic  which  is  estimated  by  the  U.S.  Geological  Survey  (USGS)  to  hold 

as  much  as  one-quarter  of  the  world’s  remaining  undiscovered  oil  and  gas  deposits 

(Borgerson  2008).  This  claim  was  rejected  by  the  UN  in  2001  based  on  a  lack  of 

evidence  provided,  which  the  Russians  have  been  working  towards  ever  since.  The 
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Canadian  government  is  particularly  opposed  to  the  Russian  claim  based  on  their  own 
assertion  that  part  of  the  Lomonosov  Ridge  constitutes  an  extension  of  their  continental 
shelf.  While  both  countries  have  proceeded  with  research  intended  to  produce  evidence 
of  their  claims,  the  Russians  garnered  global  attention  in  2007  when  they  used  a  small 
submersible  vehicle  to  plant  an  encased  national  flag  along  the  ridge  at  the  North  Pole 
(Vasilyeva  2010).  Although  this  act  was  seen  as  primarily  symbolic,  it  highlights  the 
potential  for  brinkmanship  that  exists  when  nations  debate  geographic  borders. 

Additional  disputes  have  been  brought  to  the  forefront  as  sea  ice  disappears  from 
the  Arctic  and  new  pathways  open  in  this  resource-rich  environment.  The  U.S.  and 
Canada  dispute  their  delineated  boundary  in  the  Beaufort  Sea,  Denmark  and  Canada 
disagree  over  the  territorial  rights  to  Hans  Island  between  Greenland  and  Ellesmere 
Island  (Canadian),  and  Norway  and  Russia  have  long  disputed  an  area  in  the  Barents  Sea 
thought  to  contain  vast  petroleum  deposits  under  the  sea  floor  (O’Rourke  2011).  Some 
agreement  has  been  achieved  between  the  U.S.  and  Russia  over  a  contested  area  of  the 
Bering  Sea,  and  between  Norway  and  Russia  in  the  Barents  Sea,  but  international 
discourse  is  continually  contested  as  diminishing  natural  resources  from  within  the 
continental  margins  become  increasingly  more  accessible  in  the  Arctic  Ocean 
(Isachenko  v  2010). 

The  Northern  Sea  Route  (NSR)  to  the  north  of  the  Eurasian  continent  would  open 
trade  between  northern  Europe  and  northeast  Asia  at  a  potential  cost  savings  in  time  and 
money  of  almost  40%  when  compared  to  travel  through  the  Suez  Canal  (Borgerson 
2008).  This  passage  has  been  dominated  by  Russian  shipping  traffic,  likely  due  to  their 
experience  in  these  waters  and  proprietary  rights  to  mapping  documents.  The  Northwest 
Passage  (NWP),  which  could  potentially  run  through  different  routes  in  the  Canadian 
Archipelago,  has  potential  for  cost  savings  for  shipping  between  northeast  North  America 
and  northeast  Asia  as  well  as  routing  between  northern  Europe  and  northwest  North 
America  when  compared  to  canal  routes  (Suez  or  Panama)  or  cape  routes  (Good  Hope  or 
Horn).  This  path,  while  more  modest  in  its  potential  cost  savings  (~20%)  for  the 
shipping  community,  has  been  contested  by  the  Canadians  as  an  inland  waterway  in 
which  their  government  would  maintain  sovereignty  and  control  (Borgerson  2008).  The 
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U.S.  and  the  European  Union  have  led  an  opposing  view  that  the  passage  constitutes  an 
international  strait  in  which  more  freedom  of  passage  would  prevail  (O’Rourke  2011). 
While  diplomatic  progress  and  understandings  have  prevailed  and  been  sufficient  in  the 
past,  these  passages  are  becoming  more  viable  options  in  the  summer  months,  which  will 
require  organizations  such  as  the  UN  to  adjudicate  some  of  these  disputes  in  order  to 
keep  the  relative  peace  among  Arctic  Nations. 

C.  U.S.  NAVY  ARCTIC  CAPABILITIES 

The  2010  Quadrennial  Defense  Review  (QDR)  identified  climate  change  as  one 
of  several  key  geopolitical  trends  that  may  influence  future  conflict  (Office  of  the 
Secretary  of  Defense  2010).  It  also  acknowledged  the  effect  that  climate  change  would 
have  on  shaping  future  operational  environments  and  defining  new  roles  and  mission 
requirements  for  the  Department  of  Defense  (DoD).  The  impact  of  these  changes  has 
forced  the  Navy  to  assess  its  capabilities  and  identify  shortcomings  as  a  means  of 
adapting  to  new  battlefields  and  areas  of  interest  in  a  new  climate  regime.  To  implement 
this  new  direction,  the  Navy  published  the  Climate  Change  Roadmap  in  2010  to  outline 
current  and  desired  capabilities  in  the  realm  of  observing,  predicting,  and  adapting  to 
climate  change  issues  (TFCC  2010).  The  focus  of  this  effort,  as  an  expansion  of  the 
Arctic  Roadmap,  was  to  develop  a  five-year  plan  of  action  items  to  help  the  Navy 
determine,  through  a  capabilities  based  assessment  (CBA)  process,  how  prepared  it  was 
for  the  operational  challenges  ahead  and  where  investments  could  be  most  effectively 
made  to  prepare  for  those  challenges. 

Key  findings  from  this  effort  were  published  in  both  the  DoD  Report  to  Congress 
on  Arctic  Operations  (2011)  and  the  Arctic  Environmental  Assessment  and  Outlook 
Report  (2011).  These  reports  identified  the  operational  challenges  that  exist  in  this  harsh 
environment.  Specifically,  communications  in  the  Arctic  are  significantly  limited  due  to 
poor  satellite  connectivity,  marginal  Global  Positioning  System  coverage,  and  unreliable 
long  range  high-frequency  communication  capabilities.  These  limitations  are  caused  by 
satellite  geometry  and  magnetic  and  solar  phenomena  above  70°N  that  degrade  the 
performance  of  current  communications  systems.  Surface  vessels  proceed  with  difficulty 
in  the  Arctic  as  the  equipment  they  depend  on  can  fail  due  to  icing,  while  drifting  sea  ice 
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and  heavy  fog  conditions  impede  navigation.  Naval  air  operations  are  susceptible  to 
similar  challenges  as  ice  and  weather  conditions  limit  ship  maneuverability  and  increase 
the  difficulty  of  the  launch  and  recovery  of  aircraft  on  icy  deck  surfaces.  In  the 
submarine  environment,  aside  from  the  avoidance  of  deep  draft  sea  ice  in  regions  of 
rafting  ice  and  challenges  surfacing  in  ice  covered  regions,  the  acoustic  signature  can 
change  rapidly  as  surface  water  temperature  and  salinity  profiles  are  altered  through  the 
freeze  and  melt  cycle.  This  process  makes  acoustic  forecast  models  difficult  to  predict 
and  warming  surface  waters  can  greatly  reduce  the  range  of  detection  capabilities  (TFCC 
2011). 

These  environmental  factors  reduce  standard  operational  capabilities  for  a  Navy 
that  has  been  built  to  operate  in  much  wanner  climates.  Additional  significant  challenges 
may  arise  from  a  lack  of  infrastructure  and  planning  for  potential  operations  in  this 
environment.  Currently  there  are  no  deep  draft  ports  that  could  support  most  naval 
vessels  within  the  Arctic  basin.  The  closest  deep  water  ports  are  near  Baffin  Bay  at  the 
Air  Force  base  in  Thule,  Greenland,  and  in  the  Aleutian  Islands  at  Dutch  Harbor,  Alaska 
(Figure  8).  Land-based  support  facilities  are  also  limited  to  various  service  bases  in 
Alaska  and  Greenland,  which  leads  to  dependence  on  support  from  our  partnership  with 
Canada  for  use  of  their  bases.  In  order  to  achieve  the  goals  of  Maritime  Domain 
Awareness  (MDA)  in  the  Arctic,  the  U.S.  depends  heavily  on  remote  observations  from 
satellites  and  sporadic  data  from  submarines  and  aircraft  overflight  of  the  region  in  order 
to  maintain  observational  coverage  for  its  mission  of  Intelligence,  Surveillance,  and 
Reconnaissance  (ISR).  Air  and  Missile  Defense  (AMD)  and  Ballistic  Missile  Defense 
are  well  supported  missions  from  an  extensive  network  of  early  warning  radars  and 
communications  stations  in  the  region  that  are  monitored  from  facilities  in  Alaska  and 
Greenland.  This  network,  along  with  our  relationship  with  the  Canadians  as  part  of  the 
North  American  Aerospace  Defense  Command  (NORAD)  allows  the  DoD  to  effectively 
protect  against  existing  threats  in  this  realm. 
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Existing  DoD  Bases  and  Facilities  in  Alaska  and  the  Arctic 


I 


Figure  8.  Existing  DoD  Bases  and  Facilities  in  Alaska  and  the  Arctic 
(From  DoD,  2011) 

A  major  point  of  concern  for  the  Navy’s  ability  to  operate  in  the  Arctic  as  the  sea 

ice  continues  to  recede  is  the  fact  that  surface  ships  are  not  ice-strengthened.  In  this 

configuration,  Navy  ships  are  not  safely  deployable  into  first  year  ice  or  even  the 

marginal  ice  zone  without  the  support  of  an  icebreaker.  With  no  icebreakers  in  the 

inventory,  USN  surface  ships  are  dependent  upon  the  dwindling  U.S.  Coast  Guard 

(USCG)  icebreaking  capabilities.  The  Coast  Guard  maintains  two  heavy  polar 

icebreakers,  the  Polar  Star  and  Polar  Sea,  and  one  medium  polar  icebreaker,  the  Healy, 

mostly  used  for  research.  The  heavy  icebreakers  have  both  exceeded  their  intended 
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30-year  service  lives,  and  the  Coast  Guard  is  in  the  process  of  decommissioning  the  Polar 
Sea  (O’Rourke  2011).  This  reality  spreads  the  USN  and  USCG  very  thin  in  terms  of  its 
ability  to  respond  quickly  to  Search  and  Rescue  and  Humanitarian  Assistance/Disaster 
Response  (HA/DR)  missions  that  are  likely  to  increase  in  the  Arctic  in  the  coming  years. 
While  requests  have  been  made  to  Congress  to  augment  this  fleet  by  adding  up  to  six 
heavy  and  four  medium  icebreakers  to  maintain  the  continuous  presence  requirements  of 
the  Naval  Operations  Concept,  approval  has  so  far  been  denied  which  places  any 
augmentation  out  at  least  ten  years  (O’Rourke  2011). 

Finally,  a  key  challenge  identified  from  the  CBA  was  in  the  current  ability  to 
observe  and  forecast  both  weather  and  sea  ice  for  operations  in  the  Arctic.  The  Naval  Ice 
Center  in  Suitland,  Maryland  provides  sea  ice  analysis  and  forecasting  products  primarily 
for  the  submarine  community  when  they  operate  in  the  region,  but  has  the  ability  to 
support  any  mission  with  operational  analysis,  as  well  (DoD  2011).  Shortcomings  exist 
in  the  development  of  these  products,  partially  due  to  the  communication  issues  to  this 
remote  region,  but  also  due  to  a  lack  of  data  coverage  and  in  depth  understanding  of  sea 
ice  dynamics,  which  are  integral  to  the  ability  to  predict  sea  ice  motion.  While  the  U.S. 
has  made  significant  investments  into  the  collection  of  environmental  data  in  the  Arctic 
through  the  establishment  of  an  Arctic  circumpolar  observing  network  from  partnerships 
among  various  U.S.  agencies,  academic  collaborators,  and  Arctic  residents,  sea  ice 
dynamical  processes  have  been  more  difficult  to  observe  (O’Rourke  2011)  than  ice  extent 
and  motion  fields.  Satellite  imagery,  buoy  networks,  and  in  situ  observations  made  from 
research  cruises  and  ice  camps  are  able  to  provide  scientists  with  observations  of  current 
mechanics  on  the  small  regional  scales;  but  the  underlying  dynamics  need  to  be 
understood  in  order  for  numerical  modeling  to  be  able  to  accurately  predict  basin-wide 
ice  coverage  for  operational  requirements,  as  well  as  long-term  predictions. 

D.  OBSERVATIONAL  AND  TRAINING  OPPORTUNITIES 

With  a  rich  history  as  a  border  nation  to  the  Arctic  along  the  Alaskan  coast  and  a 
member  of  the  Arctic  Council,  the  U.S.  government  has  a  vested  interest  in 
environmental  changes  as  well  as  its  operational  capabilities  in  the  region.  The  recent 
pattern  of  declining  sea  ice  and  general  progression  towards  seasonal  ice  cover  in  the 
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region  has  renewed  concerns  about  territorial  claims  and  the  vast  resources  available  in 
the  area.  Understanding  the  processes  that  drive  the  amount  of  open  water  and  predicting 
the  length  of  time  that  Arctic  sea  passages  will  be  open  to  navigation  will  be  paramount 
in  determining  the  amount  of  resources  the  U.S.  Navy  will  have  to  allocate  to  that  region 
in  the  foreseeable  future. 

While  there  is  low  potential  for  anned  conflict  in  the  region  in  the  foreseeable 
future,  the  existing  defense  infrastructure  was  deemed  adequate  to  meet  current  near-  to 
mid-term  U.S.  national  security  requirements  (DoD  2011).  In  order  to  remain  prepared  to 
operate  in  the  Arctic,  the  Navy  needs  to  remain  engaged  in  current  research  and 
operational  testing  exercises  to  maintain  and  build  upon  current  capabilities  and  to 
observe  changes  in  the  Arctic  as  they  happen.  For  the  last  two  years,  the  U.S.  Navy  has 
participated  in  the  annual  Canadian  exercise,  Operation  NANOOK,  in  which  it  has  sent 
guided  missile  destroyers  to  train  for  HA/DR  scenarios  in  the  Arctic.  The  U.S.  Marine 
Corps  has  also  participated  in  the  Norwegian  invitational  exercise  COLD  RESPONSE 
and  the  Air  National  Guard  and  Army  Reserve  have  participated  in  the  Arctic  Care 
mission,  which  provides  medical  treatment  to  remove  villages  across  Alaska.  Along  with 
these  operations-based  exercises,  several  research  opportunities  can  offer  experience  to 
military  personnel  while  improving  Arctic  sea  ice  research  knowledge  within  the  DoD. 
Continued  involvement  in  the  biennial  Ice  Exercise  (ICEX)  program  run  by  Commander, 
Submarine  Force,  Pacific  offers  both  an  opportunity  for  submarine  training  of  equipment 
and  tactical  procedures,  and  provides  a  base  for  research  through  an  agreement  with  the 
Office  of  Naval  Research  (ONR).  Utilization  of  USCG  and  other  ships  of  opportunity 
research  ventures  by  ONR  science  initiatives  will  further  improve  our  knowledge  of  the 
changing  operational  environment  in  the  Arctic. 

Apart  from  its  relevance  to  the  broader  impacts  of  global  climate  change,  a  better 
understanding  of  sea  ice  dynamics  will  assist  the  Navy  in  detennining  how  to  effectively 
operate  in  such  harsh  environments.  With  the  development  of  Task  Force  Climate 
Change  and  the  Arctic  Roadmap,  the  Navy  has  positioned  itself  at  the  forefront  of 
decisions  concerning  the  national  security  interests  of  the  region.  These  concerns 
include,  but  are  not  limited  to,  safety  of  navigation,  resource  development,  vessel  design, 
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and  the  ability  to  exploit  the  environment  for  tactical  and  operational  considerations.  If 
the  action  items  from  publications  such  as  the  Environmental  Assessment  and  Outlook 
(2011)  and  the  government  accountability  office  (GAO)  report  to  Congress  (2012)  are 
adhered  to,  the  Navy  will  be  able  to  make  informed  decisions  based  on  a  synthesis  of  the 
best  scientific  information  available. 

Uncertainty  remains  about  the  rate  and  extent  of  climate  change  in  the  Arctic,  and 
the  variable  rate  of  change  may  not  unfold  in  a  linear  fashion.  The  amount  of  increased 
human  activity  in  the  region  will  attempt  to  keep  pace  with  the  changing  environment, 
which  is  why  continual  monitoring  of  sea  ice  changes  is  required  to  make  intelligent  and 
timely  investments  that  will  improve  naval  capabilities  in  the  Arctic.  In  order  to  develop 
the  partnerships  required  for  successful  integration  of  international  policy  in  the  region, 
members  of  Task  Force  Climate  Change  will  need  to  understand  the  dynamics  of  the 
region  and  the  timeline  of  the  procurement  process  that  will  determine  the  urgency  of 
action  by  policymakers. 
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III.  DESCRIPTION  OF  NUMERICAL  MODELS  AND 
OBSERVATIONAL  DATA 

A.  MODEL  DESCRIPTIONS 

In  order  to  improve  understanding  of  the  critical  physical  processes  and  feedbacks 
in  the  Arctic  that  detennine  sea  ice  variability,  numerical  model  output  has  been 
compared  to  satellite  imagery  as  a  means  of  evaluating  which  combination  of  model 
parameters  best  represents  observational  data.  This  study  is  a  comparative  analysis  of 
regional  coupled  ice  ocean  models  which  have  been  modified  with  different  atmospheric 
forcing  and  coupling  between  parameters.  These  models  were  all  run  by  the  Naval 
Postgraduate  School  Arctic  Modeling  Effort  (NAME)  and  are  described  in  detail  below. 

1.  Polar  Ice  Prediction  System  (PIPS) 

The  Polar  Ice  Prediction  System  (PIPS)  is  a  coupled  ice-ocean  model  with  a 
horizontal  grid  spacing  of  1/12°  (or  ~9km)  projected  onto  a  rotated  spherical  grid  to 
eliminate  singularity  at  the  pole  (Maslowski  et  al.  2004).  The  model  domain  incorporates 
the  Seas  of  Japan  and  Okhotsk  in  the  west,  the  sub-Arctic  North  Pacific  and  North 
Atlantic  Oceans,  the  Canadian  Arctic  Archipelago  (CAA),  and  the  Nordic  Seas  in  the  east 
(Figure  9).  This  domain  includes  all  major  oceanic  inflow  and  outflow  regions  of  the 
Arctic  Ocean,  and  all  of  the  seasonally  ice-covered  seas  in  the  Northern  Hemisphere.  An 
artificial  channel  was  introduced  into  the  domain  across  the  North  American  continent  to 
balance  the  net  northward  water  transport  through  the  Bering  Strait  (Maslowski  et  al., 
2004).  Bathymetry  north  of  64°N  is  derived  from  the  2.5km  resolution  International 
Bathymetry  Chart  of  the  Arctic  Ocean  (IBCAO)  field  interpolated  to  the  9km  grid 
(Jakobsson  et  al.  2000;  Maslowski  and  Walczowski  2002).  Bathymetry  south  of  64°N  is 
derived  from  ETOP05  at  five-minute  resolution.  The  high  horizontal  resolution  and 
large  domain  allow  for  simulations  of  realistic  exchanges  of  mass  and  heat  between  the 
Arctic  and  its  surrounding  seas  (Maslowski  et  al.  2004).  This  model  is  considered  eddy 
pennitting  as  features  down  to  ~37km  (four  grid  points)  can  be  resolved,  but  smaller 


27 


features  down  to  the  size  of  Rossby  radius  of  deformation  in  the  Arctic  of  ~  101cm 
(Maslowski  et  al.  2008)  are  not  well  represented. 
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Figure  9.  NPS  model  domain  and  bathymetry. 

The  ice  model  consists  of  a  parallel  version  of  the  Hibler-type  sea  ice  model  with 
a  viscous-plastic  rheology,  a  zero-layer  approximation  of  heat  conduction  through  ice, 
and  two  thickness  categories  (thick  ice  and  thin  ice/open  water)  (Flibler  1979;  Zhang  and 
Hibler  1997;  Maslowski  and  Lipscomb  2003).  This  ice  model  is  coupled  to  an  ocean 
model  that  consists  of  a  regional  adaptation  of  the  Los  Alamos  National  Laboratory 
(LANL)  Parallel  Ocean  Program  (POP)  model  with  a  free  surface  fonnulation  (Dukowicz 
and  Smith  1994).  The  free  surface  and  high  resolution  allow  for  use  of  an  unsmoothed 
and  realistic  bathymetry  which  was  not  possible  with  rigid  lid  models  (Dukowicz  and 
Smith  1994).  There  are  45  fixed  ocean  depth  layers  in  the  vertical  direction  with  a 
concentration  of  eight  levels  in  the  upper  50m  and  15  levels  in  the  upper  200m  to  capture 
resolution  along  the  vast  Arctic  continental  shelves  and  slopes  (Maslowski  et  al.  2004  and 
2008). 
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The  ocean  model  was  initialized  with  climatological,  three-dimensional 
temperature  and  salinity  fields  from  the  Polar  Science  Center  Hydrographic  Climatology 
(Steele  et  al.  2001)  and  integrated  for  48  years  in  a  spin-up  mode  to  produce  realistic 
ocean  circulation  at  the  time  of  integration  (Maslowski  et  al.  2008).  The  model  was  then 
integrated  with  daily  averaged  inter-annual  atmospheric  forcing  from  the  European 
Centre  for  Medium-Range  Weather  Forecasting  (ECMWF)  for  26  years  (1979  to  2004) 
reanalysis  data  (Maslowski  and  Lipscomb  2003).  This  data  set  includes  downward 
longwave  and  shortwave  radiation,  air  temperature  and  density,  specific  humidity, 
precipitation,  wind  velocity  and  surface  wind  stress.  River  run  off  is  also  included  in  the 
ocean  model  to  allow  for  adjustment  of  the  freshwater  budget  (Maslowski  and  Lipscomb 
2003).  The  sea  ice  model  was  initialized  with  a  uniform  2m  thick  layer  of  ice  in  January 
over  an  ocean  with  climatological  surface  temperatures  less  than  0°C  to  avoid  influx  of 
freshwater  from  melting  ice  in  regions  where  sea  ice  is  not  normally  present  (Maslowski 
et  al.  2004).  Further  description  of  model  parameters  has  been  provided  elsewhere 
(Maslowski  and  Lipscomb  2003;  Maslowski  et  al.  2004). 

2.  Parallel  Ocean  Program  (POP) — Community  Ice  CodE  (CICE) 

The  POP-CICE  model  mirrors  the  PIPS  ocean  model,  domain,  bathymetry,  and 
atmospheric  forcing,  but  differs  greatly  in  its  approach  to  sea  ice  modeling.  The  CICE 
ice  model  was  developed  at  LANL  in  an  effort  to  provide  a  computationally  efficient  sea 
ice  model  to  be  run  with  fully  coupled  global  climate  models  (Hunke  and  Lipscomb 
2001).  It  can  also  run  in  stand-alone  mode  or  coupled  to  regional  models,  as  has  been 
done  in  this  study.  It  was  designed  specifically  to  be  compatible  with  the  POP  ocean 
circulation  model,  and  has  several  interacting  components.  These  components  include  a 
multi-layer  thermodynamic  model  to  compute  the  local  growth  rate  of  snow  and  ice,  an 
ice  dynamics  model  which  predicts  the  velocity  field  based  on  ice  strength,  a  transport 
model  that  describes  advection  of  ice  concentration  and  volume,  and  a  ridging 
parameterization  that  transfers  ice  among  thickness  categories  based  on  energetic 
balances  and  stain  rates  (Hunke  and  Lipscomb  2006). 

The  primary  difference  from  PIPS  ice  model  is  that  each  grid  cell  in  the  CICE 

model  contains  multiple  ice  thickness  categories  and  incorporates  elastic-viscous-plastic 
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(EVP)  dynamics,  which  is  important  for  resolving  ice  defonnation  processes  (Hunke  and 
Dukowicz  1997).  This  dynamic  model  reduces  to  a  viscous  plastic  configuration  for  long 
scale  atmospheric  forcing,  but  uses  elastic  waves  on  short  time  scales  to  improve 
response  times  and  increase  computational  efficiency  (Maslowski  and  Lipscomb  2003). 
The  most  significant  improvement  comes  from  the  multiple  ice  thickness  categories 
allowed  in  this  ice  model.  The  model  includes  five  World  Meteorological  Organization 
ice  thickness  categories,  each  having  four  vertical  ice  layers  and  one  snow  layer  (Hunke 
and  Lipscomb  2001).  Ice  strength  is  one  of  the  key  parameters  that  detennine  the  amount 
of  sea  ice  deformation  that  is  possible  and  is  dependent  upon  ice  thickness  distribution 
(Maslowski  and  Lipscomb  2003). 

This  study  will  analyze  model  output  using  different  parameterizations  and 
forcing  from  within  this  POP-CICE  model,  which  will  be  described  in  the  results  section. 
Additional  details  on  this  particular  sea  ice  model  can  be  found  in  the  CICE  User’s 
Manual  (Hunke  and  Lipscomb  2006). 

3.  Regional  Arctic  Climate  System  Model  (RACM) 

The  Regional  Arctic  Climate  System  Model  (RACM)  is  a  fully-coupled  land, 
atmosphere,  sea  ice,  and  ocean  model  developed  to  address  deficiencies  in  and  provide 
guidance  to  GCMs  in  representing  climate  change  in  the  Arctic  (Maslowski  et  al.  2012). 
RACM  was  created  under  the  support  of  the  Department  of  Energy  Earth  System 
Modeling  program  and  contains  updated  versions  of  POP  ocean  and  CICE  sea  ice  models 
described  previously  (Hunke  and  Lipscomb  2008).  Where  RACM  differs  however,  is 
that  it  incorporates  a  version  of  the  National  Center  for  Atmospheric  Research  (NCAR) 
Weather  Research  and  Forecasting  (WRF)  atmospheric  model  (Skamarock  et  al.  2005) 
that  has  been  optimized  for  polar  regions  (Hines  and  Bromwich  2008;  Bromwich  et  al. 
2009;  Higgins  et  al.  2012).  In  addition,  the  Variable  Infiltration  Capacity  (VIC)  model 
(Cherkauer  et  al.  2003;  Andreadis  et  al.  2009;  Bowling  and  Lettanmaier  2010)  represents 
land  surface  processes  and  hydrology.  All  climate  system  components  within  this  model 
are  fully  coupled  so  that  changes  made  in  one  field  impact  and  are  responded  to  by  every 
other  field  through  feedback  mechanisms.  This  means  that  unlike  the  previously 
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described  models,  where  the  atmospheric  forcing  was  prescribed,  the  atmosphere  in 
RACM  reacts  to  changes  in  the  sea  ice  model  and  so  forth  for  all  interrelated  model 
components. 

The  RACM  model  domain  encompasses  the  entire  pan- Arctic  region  (Figure  10), 
which  is  extended  from  previous  models  discussed.  The  active  ocean  and  ice  domains 
remain  the  same,  but  the  atmosphere  and  land  models  are  extended  to  include  all 
terrestrial  drainage  basins  that  bring  freshwater  to  the  Arctic  (Maslowski  et  al.  2012). 
Similar  to  previous  models,  the  ocean  and  sea  ice  models  use  a  horizontal  grid  spacing  of 
9km,  while  the  atmosphere  and  land  models  use  a  horizontal  grid  spacing  of  50km. 


RACM  Pan-Arctic  domains 


120E  105E  90E  75E  60E  45E  30E 


165W  150W  135W  120W  105W  90W  75W 


- Extended  ocean  - Ocean/sea  ice  - Arete  system 

Depth  (m) 
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Figure  10.  RACM  pan- Arctic  model  domain.  WRF  and  VIC  model  domains 

include  the  entire  colored  region.  POP  and  CICE  domains  are  bound  by 
the  inner  blue  rectangle.  Arctic  System  domain  is  bound  by  the  red  line 
(From  Maslowski  et  al.  2012,  in  press). 
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The  focus  of  the  RACM  modeling  effort  has  been  to  improve  the  realism  in 
representing  critical  small-scale  processes  and  pan-Arctic  climate  dynamics,  and  to  test 
new  physical  parameterizations,  which  will  help  guide  future  observational  requirements. 
Evaluation  of  earlier  high-resolution  NPS  models  has  shown  improvements  in  the 
seasonal  cycle  of  sea  ice  edge,  thickness,  and  volume  fluctuations  (Maslowski  et  al.  2007; 
Clement  et  al.  2005;  McGeehan  and  Maslowski  2011).  Melt  rates  from  these  earlier 
results  (Maslowski  et  al.  2007)  compare  well  with  recent  observational  data  (Kwok  and 
Cunningham,  2008;  Kwok  et  al.,  2009),  which  suggest  that  near  ice-free  summers  in  the 
Artie  may  occur  in  the  near  future  (Maslowski  et  al.  2012).  These  results  provide 
motivation  to  utilize  the  advanced  coupling  mechanism  of  a  high-resolution  RACM 
model  to  more  accurately  reproduce  the  sea  ice  dynamics  that  influence  future  climate 
predictions. 

This  thesis  further  analyzes  the  synthesis  of  the  coupled  RACM  model  output  as  a 
means  of  comparison  to  both  prior  modeling  efforts,  and  observational  data  from  satellite 
imagery.  While  the  coupling  of  constituents  found  in  the  RACM  model  does  not  offer 
direct  comparison  to  models  with  prescribed  realistic  atmospheric  forcing  (PIPS  and 
POP-CICE),  the  qualitative  comparison  to  observational  data  is  appropriate  for  further 
model  validation.  Additional  details  of  the  RACM  model  are  provided  elsewhere 
(Higgins  et  al.  2012;  Maslowski  et  al.  2012). 

B.  RADARSAT  GEOPHYSICAL  PROCESSOR  SYSTEM  (RGPS)  DATA  SET 

RADARSAT  Geophysical  Processor  System  (RGPS)  is  a  computer  system  that 
analyzes  images  of  sea  ice  in  the  Arctic  Ocean  obtained  from  RADARSAT- 1  as  part  of 
NASA’s  MEaSUREs  (Making  Earth  Science  data  records  for  Use  in  Research 
Environments)  program.  This  MEaSUREs  project  has  sought  to  compile  data  across 
several  different  Earth  Science  fields  to  produce  a  unified  and  coherent  data  set  to  address 
specific  scientific  questions  (Jet  Propulsion  Laboratory  (JPL)  2012).  RGPS  data  is 
obtained  by  tracking  sea  ice  motion  from  satellite  images  on  a  high-resolution  grid,  which 
is  developed  into  data  sets  of  small-scale  kinematic  and  deformation  fields.  These  fields 
have  been  obtained  as  a  primary  means  of  obtaining  basin-wide  observations,  improving 
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model  sea  ice  dynamics  (Coon  et  al.  2007),  and  to  help  detennine  the  impact  of  Arctic 
sea  ice  distribution  and  deformation  on  overall  climate  stability  (Lindsay  and  Stern 
2003). 

RADARSAT-1  was  developed  by  the  Canadian  Space  Agency  (CSA)  and 
launched  by  the  National  Aeronautics  and  Space  Administration  (NASA)  in  November 
1995  in  exchange  for  access  to  the  data  produced  by  the  satellite  (CSA  2012).  The 
satellite  carries  onboard  a  synthetic  aperture  radar  (SAR)  that  is  focused  on  the  Arctic 
region  for  specific  tasking  from  the  RGPS  project  to  detect  sea  ice  motion  and 
deformation  processes  (JPL  2012).  The  SAR  sensor  transmits  microwave  energy  at  5.3 
GHz  in  the  C  band  wavelength  (5.6cm).  The  radar  returns  are  recorded  and  processed 
into  imagery  at  the  satellite  facilities  in  Alaska  and  California  (CSA  2011).  Using  this 
active  imaging  technique  allows  mapping  to  be  produced  for  analysis  both  day  and  night 
as  well  as  through  various  atmospheric  conditions  common  to  the  region  such  as  cloud 
cover  and  snow  storms  (CSA  2012). 

Coverage  of  the  Arctic  occurs  as  the  satellite  orbits  the  Earth  ~  14  times  a  day  with 
a  period  of  100.7  minutes.  A  retrograde  inclination  (98.6  degrees)  determines  the  path 
around  the  Earth  in  opposite  directional  rotation  as  the  satellite  precesses  one  degree  per 
day  to  keep  up  with  the  Earth’s  rotation  around  the  Sun.  This  allows  for  an  exact  repeat 
path  flown  every  24  days  and  complete  coverage  of  the  Arctic  approximately  every  6 
days  (CSA  2010).  Since  RADARSAT-1  is  a  right-looking  satellite,  it  ascends  from  the 
equator  facing  east,  and  descends  from  the  North  Pole  facing  west  which  enables  images 
to  be  taken  from  opposite  sides  of  the  features  on  the  surface.  This,  along  with  the  ability 
to  alter  the  direction  and  width  of  the  beam  path,  provides  observation  from  many 
different  perspectives  as  the  satellite  tracks  different  features  (CSA  2012). 

As  with  all  non-polar  orbiting  satellites,  RADARSAT-1  leaves  a  hole  in  coverage 
over  the  central  Arctic  since  the  orbit  and  range  of  the  SAR  never  quite  reaches  the  pole 
(Figure  11).  This  gap  is  relatively  small  in  comparison  to  other  satellite  constellations, 
however,  as  the  areal  extent  missed  by  RGPS  is  only  about  a  two-degree  radius  (~200km) 
in  wide  beam  mode  (Lindsay  and  Stern  2003).  Lack  of  coverage  in  this  area  does  not 
greatly  impact  deformation  analysis  of  sea  ice  as  predominant  ice  motion  features  occur 
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throughout  the  Arctic  basin  and  are  often  concentrated  in  regions  with  thinner  ice 
coverage  than  is  found  at  the  pole.  With  an  overall  decrease  in  ice  thickness  and 
advection  of  deformations  into  this  region,  extrapolation  of  nearby  features  is  required  to 
cover  this  gap  in  the  data  for  comparison  with  complete  Arctic  coverage  found  in  model 
data  (Lindsay  and  Stern  2003).  For  the  purposes  of  this  research  however,  the  region 
around  the  pole  is  avoided  for  direct  comparison  between  RGPS  and  model  results. 


Figure  1 1 .  Nominal  repeat  coverage  of  the  Arctic  Ocean  over  a  24-day  cycle 
(From  Kwok  and  Cunningham  2000) 

In  wide  beam  mode,  the  SAR  provides  a  swath  of  observations  460  kilometers 
wide  with  a  pixel  resolution  of  100  meters  (Lindsay  and  Stern  2003).  This  allows  for 
reasonable  acuity  to  distinguish  different  sea  ice  deformation  parameters  within  the 
covered  area.  At  high  latitudes  in  the  Arctic  there  is  substantial  overlap  of  swath  paths, 
which  can  be  beneficial  to  tracking  of  distinct  features,  but  has  to  be  filtered  out  during 
analysis  so  the  same  feature  is  not  accounted  for  multiple  times  which  would  bias  the 
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data  significantly  (Lindsay  and  Stem  2003).  Under  RGPS  guidance,  complete  coverage 
of  the  Arctic  is  scheduled  for  every  three  days  and  this  is  generally  attained  in  the  western 
Arctic  using  the  Alaskan  station  mask,  which  covers  the  Beaufort,  Chukchi  and  East 
Siberian  Seas.  The  Tromso  mask  covers  the  eastern  Arctic,  including  the  Laptev,  Kara, 
and  Barents  seas  (including  Fram  Strait)  and  is  generally  imaged  every  six  days  due  to 
different  tasking  in  that  region  (Lindsay  and  Stern  2003).  A  sampling  period  of  less  than 
three  days  would  be  more  optimal  to  capture  small-scale  kinematic  features  that  occur 
within  that  interval  (Kwok  2006),  but  three-day  coverage  in  the  western  Arctic  has 
provided  for  observational  analysis  of  ice  motion  that  was  not  possible  before  the  launch 
of  RADARS  AT- 1 . 

The  data  set  collected  from  RGPS  represents  one  of  the  most  comprehensive  sets 
of  observations  in  the  harsh  environment  that  is  the  Arctic.  The  data  is  read  and  analyzed 
at  JPL  in  Pasadena,  California,  as  well  as  at  the  Alaska  SAR  Facility  (ASF)  in  Fairbanks, 
Alaska  (Kwok  and  Cunningham  2000).  This  redundancy  has  allowed  for  both  validation 
of  the  data  and  provided  assurance  that  data  is  received  from  the  satellite  in  the  case  that 
either  site  is  down  for  any  reason.  Understanding  the  accuracy  of  the  data  is  imperative 
to  the  ability  to  justify  the  analysis  and  any  conclusions  drawn  from  its  research. 

1.  RGPS  Data  Validation 

While  RGPS  data  has  been  developed  into  several  products  that  address  sea  ice 
kinematics,  validation  and  analysis  for  this  study  will  be  focused  on  the  more 
fundamental  products  of  Lagrangian  ice  motion  (LIM)  and  sea  ice  deformation  fields. 
The  RGPS  data  set  is  developed  with  a  Lagrangian  approach  in  which  sea  ice  velocity 
fields  are  calculated  from  trajectories  of  tracked  ice  features  from  an  initial  unifonnly 
spaced  grid  of  ~10  km  (Kwok  et  al.  2008).  The  initial  grid  then  moves  and  deforms 
throughout  the  selected  time  frame  to  produce  displacement  and  trajectory  fields  of  the 
ice  pack.  An  ice  motion  tracking  algorithm  was  developed  at  ASF  to  achieve  the  ability 
within  the  RGPS  program  to  pair  images  from  successive  satellite  passes  and  match 
features  within  these  images  (Kwok  et  al.  1990).  The  tracking  device  produces  an  array 
of  points  from  the  center  of  a  grid  point  of  interest  at  an  initial  time,  and  acquires  the 

trajectory  for  each  point  as  it  moves  with  the  ice  field  (Kwok  et  al.  1998).  This  allows  for 
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tracking  of  ice  floes  to  produce  velocity  fields  from  which  deformation  products  such  as 
divergence,  vorticity,  and  shear  can  be  derived. 

RGPS  data  has  been  validated  through  a  number  of  methods.  Initially, 
comparisons  were  made  by  analyzing  the  velocity  fields  produced  by  the  two  processing 
sites  at  JPL  and  ASF.  Through  use  of  the  same  tracking  algorithms,  automation  of  the 
velocity  fields  by  both  sites  has  produced  minimal  error  in  which  the  average 
displacement  difference  was  90  meters,  which  is  smaller  than  the  pixel  size  of  the  images 
(Lindsay  and  Stern  2003).  More  significant  error  is  introduced  however  when  there  is 
ambiguity  in  the  equivalent  backscatter  within  the  images  produced  (Kwok  and 
Cunningham  2002).  An  automatic  feature  within  the  system  flags  areas  in  which  there  is 
low  contrast  between  successive  images  or  excessive  defonnation  and  rotation  over  the 
variable  time  step.  In  this  scenario,  operator  intervention  is  required  to  determine  the 
discrepancy  and  make  a  judgment  on  the  outcome.  From  this  personal  interpretation, 
errors  are  introduced  on  the  order  of  ~2.5  km  from  initial  position  of  a  feature  on  the  first 
image  to  the  determined  position  on  the  second  image  (Lindsay  and  Stem  2003).  Manual 
interpretation  is  only  required  however  in  less  than  0.5%  of  the  cases  analyzed  for 
validation.  Therefore,  from  this  perspective  the  error  is  assumed  to  be  no  more  than  the 
pixel  size  of  the  images,  or  100  m,  which  is  consistent  with  the  observed  95th  percentile 
of  the  displacement  differences  (Lindsay  and  Stern  2003). 

RGPS  data  error  has  also  been  validated  through  comparison  with  the  more 
limited  data  set  from  buoy  trajectories  from  the  International  Arctic  Buoy  Program 
(IABP).  This  program  maintains  a  network  of  buoys  in  the  Arctic  to  provide  real  time 
operational  and  research  data  to  interested  users.  The  position  error  for  these  buoys  has 
been  estimated  as  0.3  km  by  Thorndike  and  Colony  (1980).  Correlation  of  the 
displacement  vectors  between  buoy  and  RGPS  data  was  performed  by  Lindsay  and  Stem 
(2003)  and  determined  to  be  0.996  for  all  cases  and  0.998  in  cases  without  manual 
intervention.  RGPS  locational  uncertainty,  when  compared  to  in  situ  buoy  position  error, 
can  therefore  be  assumed  to  be  on  the  same  order  as  buoy  estimates,  or  approximately 
323  meters  (Lindsay  and  Stern  2003). 


36 


The  remaining  errors  to  be  assessed  arise  from  conflict  in  geolocation  of  points 
and  those  errors  that  propagate  through  deformation  calculations.  For  deformation 
analysis,  geolocation  errors  do  not  play  a  role  in  the  statistics  since  displacements  are 
calculated  from  relative  motion  vectors  from  a  common  set  of  images,  so  these  errors  can 
be  ignored  (Lindsay  and  Stern  2003).  Deformation  calculations  are  derived  from  three 
invariants  detennined  from  the  spatial  gradient  of  the  velocity  field.  These  invariants  of 
divergence,  shear,  and  vorticity  do  not  change  relative  to  the  coordinate  system  chosen, 
which  allows  them  to  be  utilized  over  any  grid.  The  errors  in  calculating  these 
deformations  are  dependent  upon  tracking  errors  and  the  amount  of  time  between  images. 
Assuming  a  time  differential  of  three  days,  the  error  for  all  three  invariants  is  estimated  to 
be  0.5%  per  day  by  Lindsay  and  Stem  (2003).  This  error  in  defonnation  calculation  can 
be  reduced  by  grouping  cells  and  taking  an  areal  average,  but  this  comes  at  the  expense 
of  spatial  resolution. 

This  unprecedented  observational  coverage  has  given  scientists  a  unique  look  into 
the  temporal  changes  in  the  ice  pack,  which  will  offer  a  better  understanding  of  the 
momentum,  mass,  and  energy  balance  of  the  Arctic  system.  To  date  this  data  has  been 
used  extensively  in  scientific  studies  of  the  region  for  a  variety  of  analysis  including  sea 
ice  areal  coverage,  production,  thickness,  and  defonnation. 

2.  RGPS  Sea  Ice  Deformation  Calculations 

Sea  ice  velocity  is  the  fundamental  observational  field  used  by  RGPS  to  calculate 
deformation  characteristics  of  the  Arctic  sea  ice.  Sea  ice  motion  is  calculated  from  time- 
sequential  satellite  imagery  taken  from  RGPS  data  sets  (Kwok  2010).  By  observing  how 
the  identifiable  features  in  the  ice  pack  move  from  one  observation  to  the  next,  average 
ice  velocity  gradients  are  derived  from  cross-conelation  between  snapshots  of  the  same 
feature.  Using  a  Lagrangian  approach,  an  array  of  trajectories  is  produced  from 
movement  of  individual  ice  features  in  space  and  time  (Kwok  et  al.  1995),  from  which 
deformation  fields  can  be  derived. 

Deformation  within  sea  ice  is  derived  from  strain  rates  caused  by  spatial 
variability  in  the  ice  motion  fields  (Lindsay  and  Stern  2003).  There  are  three  invariants 

that  can  be  calculated  from  the  velocity  gradients  formed  by  moving  ice  fields.  The  first 
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is  divergence,  which  is  a  measure  of  the  rate  of  area  change  of  the  ice.  The  second  is 
vorticity,  which  is  a  measure  of  the  rotation  rate.  The  third  is  shear,  which  is  a  scalar 
measure  of  the  relative  deformation  of  adjacent  surfaces  (Kwok  2006).  These  invariants 
are  calculated  by  the  following  equations  (Lindsay  and  Stern  2003;  Kwok  et  al.  2008): 


_ .  du  ,  dv 

Diverqence  — - 1 - 

a  dx  dy 

Vorticity  = - 
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The  strain  rate  components  are  computed  by  RGPS  from  approximations  of  the  line 
integral  around  the  boundary  of  each  of  the  cells.  For  example, 


du 

dx 


j$u-dy  =  ^X?=i(ui+i  +  Ui)(yi+1-yi) 


(4) 


where  n  is  the  number  of  vertices  for  the  cell  and  A  is  the  area  of  the  call.  The  area  is 
computed  by, 


A  =  ±Z?=1(xtyt+1-yiXi+1) 


(5) 


where  x;  and  y;  are  the  position  coordinates  of  the  vertices.  In  a  similar  manner,  the  rest 
of  the  gradients  are  calculated  as  follows: 


du 

dy 
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(6),  (7),  (8) 


In  this  study,  the  invariants  of  divergence  and  shear  will  be  the  focus  of  the 

comparative  analysis  between  RGPS  and  model  output.  This  study  follows  the 

framework  of  Kwok  et  al.  (2008)  in  its  approach  to  determine  what  improvements  have 
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been  made  due  to  recent  model  development.  Seasonal  drift  and  regional  deformation  of 
sea  ice  from  RGPS  fields  will  be  compared  to  several  different  runs  of  NPS  model  output 
both  qualitatively  and  quantitatively.  Seasonal  and  regional  deformation  calculations  will 
be  made  for  comparison  with  RGPS  values  as  well.  Finally,  shear  and  divergence 
snapshots  will  be  taken  to  compare  model  capability  in  representing  the  linear  kinematic 
features  (LKFs)  seen  in  RGPS  output.  Throughout,  additional  approaches  will  be 
analyzed  to  detennine  if  alternate  calculations  of  model  output  better  represents  RGPS 
analysis  and  assumptions. 
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IV.  RESULTS 


In  this  chapter,  results  of  sea  ice  drift  and  defonnation  from  the  previously 
described  models  are  compared  with  RGPS  derived  ice  motion  and  deformation  fields.  A 
more  detailed  analysis  and  discussion  of  these  results  and  their  implications  for  model 
improvements  and  future  studies  can  be  found  in  following  chapters. 

A.  SEASONAL  DRIFT  AND  DEFORMATION 

In  order  to  assess  model  performance  for  comparison  to  the  RGPS  data  set,  spatial 
gradients  and  ice  motion  are  analyzed  in  this  section  initially  on  seasonal  and  regional 
scales  before  moving  to  comparisons  of  smaller  scale  deformation  processes.  Adapted 
from  the  approach  taken  in  the  Kwok  et  al.  (2008)  paper,  seasonal  RGPS  coverage  of 
Arctic  sea  ice  is  broken  into  five  regions  based  on  average  thickness  distribution  and 
common  forcing  influences  within  each  region  (Kwok  et  al.  2008).  The  regions  are 
defined  following  Kwok  et  al.  (2008)  as  S 1  and  S2  containing  ice  from  the  Beaufort  and 
Chukchi  Seas  respectively,  which  are  covered  with  predominantly  seasonal  and  thinner 
sea  ice.  S3  contains  a  mixture  of  ice  types  with  seasonal  ice  toward  the  edges  of  the  East 
Siberian  and  Laptev  seas.  S4  and  S5  contain  generally  older  and  thicker  perennial  ice 
from  the  central  Arctic  and  Canadian  Basin,  respectively.  S4  is  the  only  region  not 
subject  to  direct  coastal  influences,  and  the  southern  edges  of  SI,  S2,  and  S3  are 
generally  exposed  to  open  ocean  in  the  fall  (Kwok  et  al.  2008). 

1.  PIPS  and  CICE  Comparison  to  RGPS  Data 

Figures  12  through  15  show  selected  sea  ice  drift  and  regional  deformation 
analysis  over  four  winter  seasons  (1996  to  2000)  of  RGPS  observations  and  various 
model  runs  based  on  availability.  In  each  figure,  the  first  two  panels  show  the  initial 
coverage  and  boundaries  of  RGPS  observations  in  the  fall  and  the  advected  and  defonned 
regions  at  the  end  of  April  from  RGPS  ice  drift,  with  average  displacement  represented 
by  vectors  in  the  center  of  each  region  (in  color).  Following  panels  show  various  model 
runs  and  their  derived  ice  motion  fields  for  the  same  regional  boundaries  at  the  end  of 
April  for  comparison  to  RGPS  results.  Model  ice  drift  regional  boundaries  and  their 
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associated  average  displacement  vectors  are  shown  in  black  for  comparison  to  RGPS 
results.  Model  results  are  derived  from  daily  average  sea  ice  velocity  fields  from 
Eulerian  model  output.  Sampling  of  RGPS  ice  motion  estimates  are  not  uniform  in  time 
and  space,  so  model  output  is  transposed  into  Lagrangian  trajectories  to  match  RGPS 
output  (Kwok  et  al.  2008).  Additional  details  of  this  process  can  be  found  in  the  RGPS 
Data  User’s  Handbook  (Kwok  and  Cunningham  2000). 

The  PIPS  model  output  (noted  as  NPS  in  the  figure)  was  analyzed  as  in  the  Kwok 
et  al.  (2008)  paper  to  ensure  continuity  in  techniques  that  were  used  in  this  thesis  for 
deriving  motion  fields.  Several  runs  of  the  POP-CICE  model  were  analyzed,  two  of 
which  are  provided  here  based  on  availability  and  similar  time  frame  coverage.  Both 
runs  utilized  ECMWF  reanalysis  (ERA-15:  1979-1993)  extended  with  operational 
products  (after  1993)  as  atmospheric  forcing.  (An  extensive  description  of  the  ERA- 15 
can  be  found  at:  http://www.ecmwf.int/research/era/ERA- 1 5 D.  This  dataset  has  been 
used  widely  in  the  modeling  community  as  a  consistent  dataset  for  various  hindcast 
studies.  The  POP-CICE  R2  (noted  as  NPS  CICE  R2)  model  run  modified  the  sea  ice 
strength  parameterization  based  on  the  relationship  of  ice  thickness  distribution  and  the 
behavior  of  sea  ice  (Rothrock  1975;  Hibler  1980;  Hopkins  and  Hibler  1991;  Zhang  and 
Rothrock  2003;  Maslowski  and  Lipscomb  2003).  The  POP-CICE  S3  (noted  as 
NPS  CICE  S3)  model  run  utilized  the  same  ice  strength  parameterization,  but  also 
altered  coupling  parameterizations  for  calculating  stresses  from  the  atmosphere  and  ocean 
to  the  sea  ice  (Hunke  and  Lipscomb  2006). 
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Figure  12.  Differences  in  regional  advection  and  deformation  between  model  and  RGPS  estimates,  November  1996  to  April  1997. 
Model-derived  regional  boundaries  are  in  black,  RGPS  boundaries  are  in  color.  Vectors  (model  in  black,  RGPS  in 
color)  near  the  center  of  each  region  show  the  average  displacement  of  the  boundary  samples  computed  from  model 
and  RGPS  data. 
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Figure  13.  Differences  in  regional  advection  and  deformation  between  model  and  RGPS  estimates,  November  1997  to  April  1998. 
Model-derived  regional  boundaries  are  in  black,  RGPS  boundaries  are  in  color.  Vectors  (model  in  black,  RGPS  in 
color)  near  the  center  of  each  region  show  the  average  displacement  of  the  boundary  samples  computed  from  model 
and  RGPS  data. 
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Figure  14.  Differences  in  regional  advection  and  deformation  between  model  and  RGPS  estimates,  October  1998  to  April  1999. 

Model-derived  regional  boundaries  are  in  black,  RGPS  boundaries  are  in  color.  Vectors  (model  in  black,  RGPS  in 
color)  near  the  center  of  each  region  show  the  average  displacement  of  the  boundary  samples  computed  from  model 
and  RGPS  data. 
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Figure  15.  Differences  in  regional  advection  and  deformation  between  model  and  RGPS  estimates,  November  1999  to  April  2000. 
Model-derived  regional  boundaries  are  in  black,  RGPS  boundaries  are  in  color.  Vectors  (model  in  black,  RGPS  in 
color)  near  the  center  of  each  region  show  the  average  displacement  of  the  boundary  samples  computed  from  model 
and  RGPS  data. 
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The  displacement  and  deformation  of  these  regional  sea  ice  boundaries  in  the 
Arctic  represent  the  large-scale  drift  patterns  and  response  of  sea  ice  to  the  prevailing 
forcing  mechanisms  over  a  winter  season.  These  results  have  been  quantified  and 
summarized  in  the  following  table  to  help  determine  the  correlation  between  RGPS  and 
model  results,  as  well  as  to  determine  if  different  parameterizations  within  individual 
model  runs  have  made  improvements  to  previous  modeling  efforts.  Table  1  displays 
calculations  of  the  total  displacement  of  individual  regions  from  the  different  model  runs 
as  a  fraction  of  RGPS  displacement,  the  directional  difference  between  model  and  RGPS 
displacement  vectors,  and  the  percentage  of  regional  area  change  of  model  output 
compared  to  RGPS  estimates  (Table  1). 
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Table  1 .  Average  seasonal  RGPS  displacements  (Drgps)  of  the  boundary  points  of  the  five  regions,  model  displacements 

(Dmodel)  as  a  fraction  of  Drgps,  directional  differences  (A9)  between  model  and  RGPS  displacement  vectors  (negative 
to  the  right  of  RGPS  vector),  seasonal  percentage  area  change  of  the  live  regions,  and  model  area  change  difference  from 
RGPS  estimates. 
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2.  RACM  Comparison  to  RGPS  Data 

Using  the  same  analysis  techniques  as  in  the  previous  section,  results  from  the 
RACM  model  were  compared  with  RGPS  regional  drift  and  deformation  estimates  as  a 
means  of  determining  improvements  made  in  the  fully  coupled  model.  Since  the 
atmosphere  is  fully  predicted  in  this  model,  this  analysis  is  simply  a  statistical 
comparison  due  to  the  difference  in  atmospheric  forcing  and  sea  ice  response  found  in  the 
RACM  model.  At  the  time  of  this  publication,  RACM  model  output  is  only  available 
through  the  mid-1990s,  so  a  generic  RGPS  year  was  chosen  for  comparison  since 
observational  data  was  not  available  for  the  same  years.  The  winter  season  of  2000-2001 
was  chosen  due  to  the  anticyclonic  nature  of  the  resulting  observations,  which  matched 
closely  with  the  behavior  of  sea  ice  from  the  model  output  during  the  1991  and  1992 
winters  that  were  assessed  (Figure  16). 

Four  cases  from  the  RACM  model  output  were  analyzed  in  comparison  with 
RGPS  ice  motion  estimations.  RACMa  was  used  as  a  control  case.  This  version  of  the 
model  included  a  momentum  flux  correction  in  the  atmospheric  model  (WRF).  This  flux 
correction  required  that  the  u*  value  calculated  and  used  by  WRF  be  averaged  with  the 
u*  value  obtained  from  the  coupler  in  RACM.  In  a  perfect  coupling  case,  all  u*  values 
should  derive  from  the  coupler,  since  these  values  consistent  with  those  used  by  other 
component  models  in  RACM  (POP,  CICE  and  VIC).  The  u*  flux  correction  in  WRF  was 
found  to  be  neccessary  due  to  instabilities  resulting  from  previous  land-atmosphere 
coupling  in  RACM.  With  new  land-atmosphere  coupling  procedures,  the  RACMb  case 
sought  to  correct  the  u*  flux  correction  so  that  WRF  only  used  surface  stress  consistent 
across  the  entire  coupled  framework.  Both  RACMa  and  RACMb  model  runs  used  a  20- 
minute  dynamic  time  step  in  the  sea  ice  model.  In  the  RACMg  case,  the  dynamic  time 
step  from  was  reduced  from  20  minutes  to  five  minutes,  and  the  RACMh  case  further 
shortened  this  dynamic  time  step  down  to  one  minute.  These  reductions  in  sea  ice 
timestepping  were  based  on  analysis  of  the  paper  by  Hunke  (2001),  which  examined  EVP 
behavior  in  sea  ice  dynamic  models.  The  following  equation  was  used  to  determine  the 
impact  of  altering  the  dynamic  time  step: 
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where  P  is  the  internal  ice  pressure  (N/m),  A  is  a  function  of  the  strain  rates  (s'1),  C  is  a 
tuning  parameter  constant,  T  is  the  damping  time  scale  (s),  Ax  and  A y  represent  the  grid 
cell  lengths  (m),  and  AtJ  is  the  EVP  time  step  (s).  From  this  relationship,  it  can  be  seen 
that  a  reduction  of  the  time  step  increases  the  relative  contribution  of  the  plastic 
component  of  deformation  in  sea  ice  dynamics  by  allowing  smaller  strain  rates  to  be 
accounted  for  by  given  compressive  stresses  while  maintaining  stability  in  the  model 
(Hunke  2001).  Results  from  this  analysis  are  displayed  in  Figure  16. 
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Figure  16.  Differences  in  regional  advection  and  deformation  between  model  and 
RGPS  estimates.  RGPS  results  (top  row)  from  2000-2001  season  are 
compared  to  two  seasons  ((1990-1991  and  1991-1992)  of  each  model  run; 
RACMa  (second  row),  RACMb  (third  row),  RACMg  (fourth  row),  and 
RACMh  (fifth  row).  Model-derived  regional  boundaries  are  in  black, 
RGPS  boundaries  are  in  color.  Vectors  (model  in  black,  RGPS  in  color) 
near  the  center  of  each  region  show  the  average  displacement  of  the 
boundary  samples  computed  from  model  and  RGPS  data. 
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As  in  the  previous  analysis  for  the  POP-CICE  model,  the  RACM  model  case 
results  have  been  quantified  by  direct  comparison  of  total  displacement  differences, 
directional  displacement  differences,  and  regional  area  change  comparisons  with  RGPS 
estimates  (Table  2). 
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1.26 

4.0 

1.7 

16.1 

2.03 

-3.3 

•  1.8 

12.7 

1.29 

-9.2 

3.3 

17.7 

S4 

1.75 

-20.1 

-5.0 

-8.0 

1.35 

-17.6 

1.0 

-2.1 

1.49 

-7.5 

8.7 

5.7 

1.85 

-II. 1 

-4.3 

-7.3 

S5 

4.02 

-24.8 

-16.9 

-15.9 

3.59 

-20.0 

-6.9 

£2 

3.26 

-17.6 

-8.5 

•7-5 

2.30 

-5.9 

-16.0 

-15.0 

Table  2.  Average  seasonal  RGPS  displacements  (Drgps)  of  the  boundary  points  of  the  five  regions,  model  displacements 

(Dmodel)  as  a  fraction  of  Drgps,  directional  differences  (A9)  between  model  and  RGPS  displacement  vectors  (negative 
to  the  right  of  RGPS  vector),  seasonal  percentage  area  change  of  the  five  regions,  and  model  area  change  difference  from 
RGPS  estimates. 
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B.  SMALL-SCALE  DEFORMATIONS 

Representation  of  small-scale  kinematics  of  defonnation  processes  may  assist  in 
the  ability  of  modeling  efforts  to  better  represent  larger-scale  sea  ice  motion  and  its 
contribution  to  sea  ice  thickness  distribution  and  changes  in  the  Arctic.  In  this  section, 
comparisons  are  made  between  deformation  fields  derived  from  RGPS  and  model  output. 
Fracture  patterns  in  the  ice  cover  represent  the  mechanical  response  of  sea  ice  to 
atmospheric  and  oceanic  forcing  (Kwok  et  al.  2008).  The  ability  of  model  output  to 
produce  similar  structures  to  observed  sea  ice  deformation  fields  will  help  to  determine 
the  relative  impact  of  kinematic  forcing  to  ice  production  and  thickness  distribution  in  the 
Arctic. 

Direct  comparison  of  the  spatial  invariants  derived  from  RGPS  and  model  results 
is  difficult  based  on  the  differences  in  how  these  features  are  derived  from  each  dataset. 
RGPS  invariants  are  calculated  from  ice  motion  fields  estimated  after  features  in  the  sea 
ice  are  tracked  over  a  mosaic  of  repeated  images  taken  from  various  satellite  passes.  On 
average,  the  western  Arctic  is  covered  over  a  three-day  period,  which  means  that  velocity 
fields  derived  from  repeated  swaths  cover  a  six-day  period.  Deformation  fields  are  then 
calculated  from  these  six-day  velocity  fields.  Model  sea  ice  velocity  fields,  by  contrast, 
are  available  at  each  point  on  the  Eulerian  grid  and  at  each  time  step,  i.e.,  every  48  min, 
with  the  standard  output  archived  on  the  daily  basis.  For  the  purposes  of  this  study, 
average  velocity  fields  were  taken  from  daily  snapshots  of  model  output. 

1.  Shear  Deformation  Analysis 

RGPS  deformation  fields  show  a  high  density  of  LKFs  within  the  ice  cover  which 
vary  in  location  and  orientation  pattern  (Kwok  et  al.  2008).  In  order  to  assess  model 
capability  of  representing  these  LKFs,  a  visual  comparison  of  RGPS  and  model 
deformation  fields  has  been  presented  in  this  section.  The  initial  approach  applied  to 
achieve  this  comparison  mirrors  the  approach  taken  in  the  Kwok  et  al.  (2008)  paper, 
whereby  daily  velocity  fields  from  the  model  output  are  used  to  calculate  the  deformation 
fields  over  the  same  six-day  period  as  in  the  RGPS  coverage.  This  is  accomplished  by 
taking  daily  velocity  snapshots  from  model  output,  calculating  the  velocity  gradients  over 
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a  Lagrangian  grid,  and  then  calculating  the  average  deformation  produced  by  the  model 
over  a  defined  time  period  (6  days  in  this  analysis).  PIPS  results  are  compared  in  this 
manner,  as  in  the  paper,  with  the  addition  of  new  model  results  from  the  POP-CICE  S3 
model  run  described  previously.  An  example  from  1997  from  the  Kwok  et  al.  (2008) 
paper  is  used  for  direct  comparison  (Figure  17). 


Dec  5-1 0,  1 997  Feb  26-Mar  3,  1 997 


Percent  /  Day 


Figure  17.  Examples  showing  correspondence  between  patterns  of  6-day  shear 

deformation  from  model  fields  and  RGPS  estimates.  Linear  kinematic 
features  (indicated  by  arrows)  in  the  model  seem  to  line  up  with  patterns 
from  RGPS  results  (Adapted  from  Kwok  et  al.  2008). 
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An  increase  in  the  density  of  deformation  features  over  the  PIPS  results  can  be 
easily  seen  in  the  POP-CICE  S3  figure,  likely  due  to  the  multiple  ice  thickness  categories 
represented  in  each  grid  cell  of  this  model.  However,  since  the  model  velocity  fields  are 
available  on  a  daily  basis,  deformations  can  be  diagnosed  on  a  daily  basis  as  well  to 
analyze  ice  kinematics  at  shorter  time  frames  than  RGPS  data  allow.  Figures  18  and  19 
offer  examples  of  this  analysis  from  the  same  time  period  (1997),  in  which  the  top  figure 
shows  the  six-day  averaged  model  deformation  field,  and  the  bottom  rows  show  daily 
snapshots  of  deformation  as  represented  by  the  model.  From  this  analysis,  it  can  be  seen 
that  many  features  present  in  the  daily  model  snapshots  are  not  present  in  the  six-day 
averaged  plots  as  calculated  following  the  approach  from  Kwok  et  al.  (2008).  Based  on 
these  results,  additional  methods  of  intercomparisons  of  model  ice  deformations  with 
those  derived  from  RGPS  are  pursued  in  the  next  section. 
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Figure  18.  December  1997  shear  example  with  6-day  average  model  field  shown  (top 
row)  with  comparison  of  daily  shear  snapshots  from  the  same  6-day  period 
(bottom  two  rows).  Additional  features  present  in  daily  snapshots 
compared  to  6-day  average. 
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Figure  19.  February  1997  shear  example  with  6-day  average  model  field  shown  (top 
row)  with  comparison  of  daily  shear  snapshots  from  the  same  6-day  period 
(bottom  two  rows).  Additional  features  present  in  daily  snapshots 
compared  to  6-day  average. 
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2. 


Alternative  Model  Shear  Estimations 


In  this  section,  an  alternative  method  of  calculating  shear  deformations  from 
velocity  fields  within  the  POP-CICE  S3  model  output  is  described.  Four  years  ( 1 996— 
1999)  are  analyzed  to  match  both  the  December  and  February  results  offered  by  the 
Kwok  et  al.  (2008)  paper.  This  method  was  developed  as  a  means  of  showing  an 
additional  methodology  to  represent  the  small-scale  deformations  found  in  the  model  in  a 
way  that  possibly  helps  to  minimize  the  differences  in  model-data  estimates  of  sea  ice 
deformations  due  to  different  types  of  data  used  (i.e.,  Lagrangian-type  observations 
versus  Eulerian-type  model  output)  and  time  required  for  basin  coverage. 

In  order  to  transpose  daily  average  velocity  fields  and  their  associated  gradients 
and  deformation  fields  calculated  by  the  model  onto  six-day  estimates  by  RGPS,  daily 
sea  ice  velocity  fields  are  averaged  over  the  same  period  (6  days)  covered  in  the  RGPS 
dataset  in  the  former  method.  By  averaging  the  model  velocity  over  six  days,  the 
magnitude  of  spatial  variability  of  the  sea  ice  velocity  field  is  reduced,  which  is  reflected 
in  the  deformation  calculations  resulting  in  a  lower  density  of  deformation  features  across 
the  domain.  This  process  was  used  to  overcome  the  difficulties  of  direct  comparison 
between  RGPS  and  model  output  discussed  earlier,  and  is  one  way  of  diagnosing  the 
model  capabilities  to  reproduce  kinematic  features,  but  it  possibly  filters  out  some  of  the 
details  actually  present  in  model  output. 

To  account  for  the  loss  of  features  resolved  in  the  model  by  this  method,  average 
velocity  fields  are  taken  at  shorter  time  steps,  and  the  deformations  are  summed  instead 
of  averaged  out  which  allows  shorter-lived  features  to  be  represented.  By  summing  the 
resulting  deformations,  some  features  can  be  counted  multiple  times  leading  to 
exaggerated  defonnation  density.  Still  this  approach  allows  for  more  realistic 
representation  of  model-resolved  deformations  in  the  central  Arctic. 

Figures  20  through  27  show  results  obtained  by  taking  average  velocities  from 
model  output  at  reduced  time  interval  compared  to  the  six-day  averages  taken  in  the 
original  RGPS  comparisons.  Five-day  averages  represent  taking  the  average  model 
velocities  from  days  1-5  and  days  2-6  during  the  same  six-day  window,  and  using  those 
velocities  to  calculate  the  defonnations.  In  a  similar  manner,  four-day  averages  take  days 
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1-4,  2-5,  3-6;  three-day  averages  take  days  1-3,  2-4,  3-5,  4-6;  and  two-day  averages 
take  days  1-2,  2-3,  3-4,  4-5,  5-6.  In  order  to  detennine  an  appropriate  summation 
range,  these  time  step  averages  were  taken  and  defonnations  were  calculated  in  the  same 
manner  by  dividing  by  the  number  of  time  steps  (top  row  in  Figures  20  through  27),  and 
also  by  simply  taking  a  sum  of  the  deformations  (bottom  row  in  Figures  20  through  27) 
calculated  from  the  applicable  time  step.  The  density  of  deformations  increases,  as 
expected,  due  to  the  reduction  of  velocity  temporal  smoothing.  Outcome  of  the 
summation  of  these  results  is  described  next. 
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Figure  20.  December  1996  example  of  alternative  shear  calculations  from  POP-CICE  S3  model  output.  Top  rows  show  shear 

deformations  from  multi-day  averaged  velocity  fields,  averaged  through  division  of  the  number  of  time  steps  in  the  6- 
day  window.  Bottom  rows  show  a  summation  of  the  shear  deformations  from  multi-day  averaged  velocity  pairs  at 
different  time  steps. 
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Percent  /  Day 


Figure  21.  December  1997  example  of  alternative  shear  calculations  from  POP-CICE  S3  model  output.  Top  rows  show  shear 

deformations  from  multi-day  averaged  velocity  fields,  averaged  through  division  of  the  number  of  time  steps  in  the  6- 
day  window.  Bottom  rows  show  a  summation  of  the  shear  deformations  from  multi-day  averaged  velocity  pairs  at 
different  time  steps. 
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Percent  Day 
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Figure  22.  December  1998  example  of  alternative  shear  calculations  from  POP-CICE  S3  model  output.  Top  rows  show  shear 

deformations  from  multi-day  averaged  velocity  fields,  averaged  through  division  of  the  number  of  time  steps  in  the  6- 
day  window.  Bottom  rows  show  a  summation  of  the  shear  deformations  from  multi-day  averaged  velocity  pairs  at 
different  time  steps. 


63 


Shear 
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Figure  23.  December  1999  example  of  alternative  shear  calculations  from  POP-CICE  S3  model  output.  Top  rows  show  shear 

deformations  from  multi-day  averaged  velocity  fields,  averaged  through  division  of  the  number  of  time  steps  in  the  6- 
day  window.  Bottom  rows  show  a  summation  of  the  shear  deformations  from  multi-day  averaged  velocity  pairs  at 
different  time  steps. 
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Figure  24.  February  1996  example  of  alternative  shear  calculations  from  POP-CICE  S3  model  output. 

Top  rows  show  shear  deformations  from  multi-day  averaged  velocity  fields,  averaged  through  division  of  the  number 
of  time  steps  in  the  6-day  window.  Bottom  rows  show  a  summation  of  the  shear  deformations  from  multi-day  averaged 
velocity  pairs  at  different  time  steps. 
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Figure  25.  February  1997  example  of  alternative  shear  calculations  from  POP-CICE  S3  model  output.  Top  rows  show  shear 

deformations  from  multi-day  averaged  velocity  fields,  averaged  through  division  of  the  number  of  time  steps  in  the  6- 
day  window.  Bottom  rows  show  a  summation  of  the  shear  deformations  from  multi-day  averaged  velocity  pairs  at 
different  time  steps. 
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Figure  26.  February  1998  example  of  alternative  shear  calculations  from  POP-CICE  S3  model  output.  Top  rows  show  shear 

deformations  from  multi-day  averaged  velocity  fields,  averaged  through  division  of  the  number  of  time  steps  in  the  6- 
day  window.  Bottom  rows  show  a  summation  of  the  shear  deformations  from  multi-day  averaged  velocity  pairs  at 
different  time  steps. 
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Figure  27.  February  1999  example  of  alternative  shear  calculations  from  POP-CICE  S3  model  output.  Top  rows  show  shear 

deformations  from  multi-day  averaged  velocity  fields,  averaged  through  division  of  the  number  of  time  steps  in  the  6- 
day  window.  Bottom  rows  show  a  summation  of  the  shear  deformations  from  multi-day  averaged  velocity  pairs  at 
different  time  steps. 


68 


C.  REGIONAL  DEFORMATION  COMPARISON 

This  section  provides  regional,  small-scale  deformation  results  from  PIPS  (NPS  in 
figures)  and  POP-CICE  S3  model  runs  compared  to  defonnation  calculated  from  RGPS 
ice  motion  fields  over  three  winter  seasons  (Figures  28  through  30).  Time  series  analysis 
is  presented  from  six-day  shear  calculations  from  model-  and  RGPS-derived  velocity 
fields  separated  into  the  five  regions  described  in  section  A.  Calculations  for  this  section 
were  performed  using  the  original  six-day  averaging  approach  described  by  Kwok  et  al. 
(2008).  These  time  series  show  spatial  and  temporal  correlations  between  model  and 
RGPS  deformation  calculations  over  entire  winter  seasons.  Deformation  magnitude  and 
temporal  phase  agreement  of  model  output  and  RGPS  estimates  are  analyzed  in  the  next 
chapter.  Statistical  analyses  of  these  plots  are  provided  in  Table  3. 
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Figure  28.  Time  series  (1997-1998)  of  the  mean  6-day  regional  shear  derived  from 
model  and  RGPS  ice  drift  for  regions  S 1  through  S5  (defined  previously) 
between  November  1997  and  April  1998. 
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Shear 

1998-  1999 


Figure  29.  Time  series  (1998-1999)  of  the  mean  6-day  regional  shear  derived  from 
model  and  RGPS  ice  drift  for  regions  S 1  through  S5  (defined  previously) 
between  November  1998  and  April  1999. 
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Shear 

1999-2000 


Figure  30.  Time  series  (1999-2000)  of  the  mean  6-day  regional  shear  derived  from 
model  and  RGPS  ice  drift  for  regions  S 1  through  S5  (defined  previously) 
between  November  1999  and  April  2000. 
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RGPS 

NPS 

S3 

RGPS 

NPS 

S3 

RGPS 

NPS 

S3 

1997-1998 

1998-1999 

1999-2000 

SI 

SI 

SI 

Mean 

1.18 

0.14 

0.67 

1.08 

0.11 

0.61 

0.85 

0.22 

0.69 

Std  Dev 

0.82 

0.20 

0.51 

0.78 

0.23 

0.33 

0.66 

0.25 

0.41 

Max 

3.64 

0.69 

2.62 

3.71 

1.18 

1.20 

2.58 

0.90 

1.72 

Min 

0.23 

0.00 

0.03 

0.01 

0.00 

0.09 

0.07 

0.00 

0.19 

CORREL 

0.48 

0.58 

0.29 

0.27 

0.61 

0.67 

S2 

S2 

S2 

Mean 

0.90 

0.22 

0.54 

0.79 

0.15 

0.45 

0.80 

0.28 

0.57 

Std  Dev 

0.62 

0.36 

0.39 

0.39 

0.21 

0.30 

0.52 

0.35 

0.30 

Max 

2.51 

1.65 

1.86 

2.00 

0.82 

1.23 

1.98 

1.64 

1.14 

Min 

0.10 

0.00 

0.04 

0.22 

0.00 

0.10 

0.11 

0.01 

0.03 

CORREL 

0.66 

0.75 

0.37 

0.35 

0.50 

0.58 

S3 

S3 

S3 

Mean 

0.89 

0.19 

0.53 

0.72 

0.10 

0.39 

0.76 

0.23 

0.65 

Std  Dev 

0.48 

0.35 

0.33 

0.47 

0.12 

0.24 

0.43 

0.23 

0.36 

Max 

2.26 

1.68 

1.72 

1.88 

0.42 

0.95 

2.02 

0.89 

1.70 

Min 

0.14 

0.00 

0.07 

0.19 

0.00 

0.08 

0.07 

0.02 

0.12 

CORREL 

0.65 

0.28 

0.26 

0.34 

0.63 

0.31 

S4 

S4 

S4 

Mean 

0.74 

0.13 

0.35 

0.70 

0.11 

0.34 

0.59 

0.20 

0.38 

Std  Dev 

0.39 

0.23 

0.27 

0.38 

0.20 

0.25 

0.43 

0.19 

0.29 

Max 

1.66 

0.94 

1.14 

1.60 

1.02 

1.20 

1.83 

0.66 

1.21 

Min 

0.09 

0.00 

0.03 

0.14 

0.00 

0.05 

0.11 

0.00 

0.05 

CORREL 

0.12 

0.06 

0.05 

0.28 

0.39 

0.53 

S5 

S5 

S5 

Mean 

0.73 

0.16 

0.37 

0.72 

0.14 

0.36 

0.64 

0.16 

0.40 

Std  Dev 

0.49 

0.29 

0.28 

0.52 

0.30 

0.22 

0.43 

0.27 

0.27 

Max 

2.01 

1.54 

1.30 

1.78 

1.59 

1.01 

1.98 

1.10 

1.02 

Min 

0.13 

0.00 

0.03 

0.05 

0.00 

0.09 

0.15 

0.00 

0.08 

CORREL 

-0.06 

0.46 

0.33 

0.53 

0.75 

0.76 

Table  3.  Statistical  comparison  of  seasonal  shear  values  from  Figures  28  through  30. 

Mean,  Standard  Deviation  (Std  Dev),  Maximum  (Max),  Minimum  (Min), 
and  Correlation  (CORREL)  values  of  shear  are  displayed  by  region  (all 
values  *10'2). 
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V.  DISCUSSION  AND  SYNTHESIS  OF  RESULTS 


Results  and  analyses  of  NPS  model  output  and  RGPS  datasets  were  presented  in 
the  previous  chapter,  with  the  primary  purpose  of  determining  model  skill  in  representing 
the  variability  in  sea  ice  kinematics  and  deformation  processes  that  impact  ice  thickness 
distribution  in  the  Arctic.  These  results  are  intended  to  offer  an  improved  insight  into 
potentially  critical  physical  processes  and  feedbacks  that  detennine  the  survivability  of 
sea  ice  in  the  Arctic.  The  following  discussion  and  conclusions  attempt  to  define  what 
improvements  have  been  made  through  sea  ice  model  parameterizations  and  forcings,  and 
how  the  performed  analyses  advance  the  simulation  and  prediction  of  ice  thickness 
distribution  in  the  Arctic.  The  challenges  of  direct  comparison  between  model  results  and 
RGPS  data  arise  from  the  differences  in  how  the  basin-wide  velocity  fields,  gradients, 
and  defonnations  are  derived.  The  main  difference  is  due  to  the  respective  frame  of 
reference  (i.e.,  Eulerian  for  model  output  vs.  Lagrangian  for  RGPS  data),  which  to  some 
degree  can  be  overcome  through  computational  translation  between  those  frameworks  to 
improve  spatial  and  temporal  details  available  from  both  datasets  compared  to  previous 
techniques  of  intercomparison.  The  strength  of  the  RGPS  dataset  is  that  it  allows  for 
comparative  analysis  of  modeled  small-scale  deformations  at  near  basin-scale  and  at 
synoptic,  seasonal,  and  inter-annual  time  scales.  This  allows  for  the  investigation  of 
short-  and  long-term  responses  of  modeled  sea  ice  to  applied  forcing  and  provides  a 
diagnostic  tool  to  improve  the  understanding  of  sea  ice  model  sensitivity  to  different 
parameterizations  and  coupling. 

A.  LARGE-SCALE  DISPLACEMENT  AND  DEFORMATIONS 

Results  from  the  regional  analysis  of  net  seasonal  displacement  and  deformation 
fields  of  PIPS  and  POP-CICE  model  output  (Figures  12  through  15)  shows  the  large- 
scale  response  of  sea  ice  to  prescribed  forcing.  From  the  PIPS  model  results,  it 
can  be  seen  that  seasonal  displacement  was  slower  than  RGPS  observations  and  the  areal 
deformation  of  each  region  was  subsequently  minimal  by  comparison.  This  was  likely 
due  to  the  simplified  ice  strength  parameterization  and  thermodynamics  represented 
in  the  PIPS  model.  The  POP-CICE  model  runs,  with  multi-layer  ice  strength 
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parameterizations  (R2)  and  atmosphere-ice  /  ice-ocean  coupling  parameterizations  (S3), 
show  a  large  increase  in  ice  motion,  which  was  overestimated  particularly  along  the 
Canadian  Arctic  Archipelago  (CAA)  and  Alaskan  coast.  Table  1  shows  that  the 
displacement  difference  alone  was  on  average  1.8  to  4.2  times  greater  in  the  R2  run  than 
RGPS  estimates  for  regions  with  coastal  influence  (SI,  S2,  S3,  and  S5).  This  ratio  was 
reduced  to  an  average  of  1.2  to  3.0  times  greater  in  the  S3  run  by  weakening  the  sea  ice 
coupling  with  the  atmosphere  and  the  ocean,  however,  total  sea  ice  motion  was  overall 
still  too  large.  Directional  differences  of  displacement  vectors  were  not  improved  from 
PIPS  results,  and  based  on  the  overestimated  displacement,  the  regional  deformations 
were  too  large  as  well.  These  results  point  to  high  model  sensitivity  to  ice  strength 
parameters  and  the  influence  of  modified  air-ice  and  ice-ocean  coupling.  In  the  CICE 
model,  ice  strength  of  a  grid  cell  is  a  function  of  the  fraction  of  the  thinnest  ice  category, 
which  appears  to  lead  to  overestimation  of  ice  drift  and  deformation  without  optimal 
tuning  of  sub-grid  parameterizations  controlling  ice  behavior  at  a  given  spatial  and 
temporal  resolution. 

Regional  displacement  and  deformation  fields  from  various  runs  of  the  RACM 
model  (Figure  16),  while  limited  to  statistical  comparison  due  to  the  predictive  nature  of 
atmospheric  forcing  (i.e.,  RACM  atmospheric  condition  being  different  from  the  real 
atmospheric  forcing  embedded  in  RGPS  data),  show  improvements  over  the  previous 
models.  The  control  cases  of  RACMa  and  RACMb,  with  20-minute  dynamic  time  steps 
in  the  sea  ice  model,  show  similar  overestimation  of  regional  displacements,  but  seem  to 
show  some  improvement  in  the  directional  differences  from  RGPS  output,  however  this 
could  be  a  coincident  with  the  selected  year  of  model  output.  The  RACMg  case,  which 
reduced  the  dynamic  time  step  to  five  minutes  bringing  the  sea  ice  solution  closer  to 
plasticity,  still  represents  too  much  ice  motion  particularly  along  the  CAA  (region  S5), 
where  the  average  displacement  ratio  ranges  from  1.07  to  3.26,  but  continues  the 
improvement  in  directional  changes  where  the  average  maximum  angular  difference  is 
~17  degrees  (compared  to  greater  than  20  degrees  in  previous  model  results).  Finally,  the 
RACMh  case  (one -minute  time  step)  greatly  reduces  the  displacement  ratio  down  to  a 
range  of  0.85  to  1.84  (1990-1991  case)  when  compared  to  RGPS  regional  displacements 
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(Table  2).  This  improvement  is  accompanied  by  further  improved  directional  differences 
in  the  displacement  vectors  in  which  the  average  maximum  angular  difference  is  reduced 
to  ~11  degrees,  and  the  geometrical  representation  of  the  regional  deformation  looks 
physically  similar  to  the  RGPS  estimates. 

The  high-frequency  air-ice  feedbacks  in  addition  to  spatially  high-resolution 
predicted  atmosphere  likely  play  a  role  in  the  improved  responses  seen  in  these  RACM 
model  runs.  Winds  produced  in  the  WRF  model  are  orders  of  magnitude  greater  than  the 
smoothed  winds  prescribed  by  the  atmospheric  forcing  used  in  the  POP-CICE  and  PIPS 
model  runs.  The  sea  ice  dynamic  time  step  size  allows  the  model  more  realistic 
representation  of  sea  ice  thickness  build  up  compared  to  longer  time  steps  when  the 
elastic  and  viscous  response  in  sea  ice  motion  is  overestimated  prior  to  plastic  failure. 
While  exact  replication  of  RGPS  analysis  is  not  expected,  statistical  diagnostics  of  sea  ice 
model  behavior  against  this  dataset  will  continue  to  help  modelers  better  understand  the 
sensitivity  and  impact  of  various  parameterizations. 

B.  SMALL-SCALE  DEFORMATIONS 

Linear  kinematic  features  (LKFs)  observed  in  the  RGPS  analysis  represent  the 
prevalence  of  fracture  and  failure  defonnation  that  occurs  in  sea  ice  as  blocks  of  different 
thickness  move  relative  to  one  another  based  on  the  varying  forcing  mechanisms  found  in 
the  Arctic.  Shear  defonnations  are  one  of  the  prime  examples,  which  display  this 
reaction  in  the  sea  ice  as  cracks  and  failure  lines  represent  the  dynamic  nature  of  the  ice 
sheet.  Difficulties  in  comparing  RGPS  deformations  to  model  representation  of  these 
features  are  inherent  to  the  difference  in  methodology  of  estimating  such  shear 
deformations.  As  briefly  mentioned  earlier,  RGPS  deformations  are  derived  from  repeat 
orbit  swaths  of  satellite  imagery  to  determine  displacements  which  then  are  used  to 
calculate  mean  sea  ice  drift,  velocity  gradients,  and  finally  deformations,  which  provide  a 
full  composite  coverage  of  the  western  Arctic  roughly  every  six  days.  Model  output 
provides  domain-wide  snapshots  of  ice  motion,  which  can  be  converted  into  defonnation 
fields  at  every  model  timestep.  To  minimize  the  differences  resulting  from  different 
types  of  data  (i.e.,  Lagrangian  vs.  Eulerian)  and  methods  to  estimate  basin-wide 
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deformation  fields,  alternative  techniques  have  been  evaluated  to  optimize 
intercomparison  of  deformation  features  between  the  two  data  sources. 

The  amount  of  LKFs  presented  by  RGPS  analysis  has  proven  difficult  to 
reproduce  by  the  PIPS  model,  likely  due  to  the  simplified  ice  thickness  parameterization 
used  in  the  model.  Figure  17  shows  that  the  POP-CICE  model,  run  with  the  same 
resolution  as  the  PIPS  model,  was  able  to  represent  an  increased  amount  of  shear 
deformations  and  LKFs  through  the  sub-grid  representation  of  multiple  ice  thickness 
categories  and  other  parameterizations  used  in  the  CICE  model.  Further  analyses 
(Figures  18  and  19)  however  revealed  that  the  averaging  of  model-derived  deformation 
fields  over  a  six-day  period  to  match  the  composite  coverage  of  RGPS  analysis  left  out 
many  short-lived  shear  deformations,  especially  in  the  deep  basin.  By  this  method,  only 
more  persistent  model  events  relative  to  the  six-day  averaging  remain. 

To  address  this  issue,  an  alternative  method  was  developed  in  an  attempt  to  pennit 
some  of  the  short-term  LKFs  present  in  model  output,  which  are  naturally  included  in 
RGPS  data.  This  method  reduces  the  time  period  (i.e.,  number  of  days)  over  which 
velocity  averages  are  calculated,  and  sums  the  results,  but  does  not  divide  (average)  them 
by  the  number  of  days  for  each  defonnation  increment.  A  disadvantage  of  this  approach 
is  that  at  shorter  averaging  time  periods,  persisting  defonnation  features  can  be  counted 
multiple  times  resulting  in  deformation  fields  with  regions  of  exaggerated  magnitudes. 
Figures  20  and  27  however,  show  that  shortening  the  time  step  and  summing  the 
deformation  fields  reveal  additional  features  from  the  model  output.  A  four-day 
summation  appears  to  represent  a  good  balance  between  too  many  repeat  counts  along 
the  margins,  and  higher  resolution  of  short-term  features  in  the  central  Arctic.  This  four- 
day  summation  (bottom  middle  panel  in  each  figure)  represents  velocities  averaged  from 
the  beginning  (days  1-4),  middle  (days  2-5),  and  end  (days  3-6)  of  the  six-day  period 
covered  by  RGPS  analysis.  The  areas  with  overly  increased  shear  percentages  lie 
primarily  outside  the  region  represented  in  the  RGPS  plots,  so  it  is  difficult  to  quantify 
their  realism.  While  there  is  no  optimal  comparison  of  model  output  to  RGPS  coverage, 
the  four-day  summations  display  the  full  breadth  of  features  that  are  present  in  the  model, 
and  are  in  better  agreement  with  RGPS-derived  deformation  fields. 
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C.  REGIONAL  DEFORMATION  COMPARISON 

While  it  is  not  expected  that  sea  ice  modeling  efforts  should  be  capable  of  exactly 
reproducing  events  from  the  past,  especially  at  the  synoptic  scale  of  the  atmospheric 
forcing  which  provides  most  of  the  impetus  for  sea  ice  motion,  seasonal  analysis  of 
deformation  features  can  be  expected  to  approach  the  general  timing,  if  not  the  absolute 
magnitude  of  events  based  on  observational  re-analysis.  In  this  respect,  the  POP-CICE 
S3  model  shows  a  vast  improvement  of  deformation  representation  over  the  PIPS  model 
presented  in  the  Kwok  et  al.  (2008)  paper.  Comparison  of  the  time  series  of  regional 
shear  over  the  three  seasons  presented  (Figures  28-30),  the  PIPS  model  displays  only 
occasional  periods  of  increased  shear  represented  during  large  events  in  the  RGPS  record, 
but  overall  is  significantly  underestimating  the  amount  of  shear  present  in  the  sea  ice. 
Almost  no  shear  is  present  in  PIPS  for  weeklong  to  monthly  periods  where  RGPS  shows 
an  entirely  different  behavior  of  the  sea  ice.  As  discussed  previously,  this  is  likely  due  to 
overestimation  of  ice  strength,  which  is  determined  as  a  function  of  the  mean  grid  cell 
thickness  in  PIPS  compared  to  being  a  function  of  the  fraction  of  the  thinnest  ice  in  CICE 
model.  The  POP-CICE  S3  model  results  show  an  improvement,  not  only  in  the  phase 
agreement  of  most  of  the  events,  but  also  by  approaching  the  magnitude  of  shear 
deformations  represented  in  RGPS  data. 

The  increase  in  magnitude  of  shear  events  throughout  the  seasons  analyzed  is 
expected  from  the  improved  sea  ice  model  representing  multiple  thickness  layers.  As  can 
be  seen  in  Table  3,  the  POP-CICE  S3  mean  seasonal  shear  values  for  all  regions  defined, 
range  from  two  to  seven  times  the  shear  magnitudes  represented  in  the  PIPS  model. 
While  on  average,  these  values  still  lag  behind  the  magnitude  of  shear  represented  in  the 
RGPS  data,  this  is  a  significant  improvement  especially  given  the  low-resolution  and 
daily  averaged  atmospheric  forcing  prescribed  in  those  runs.  It  is  also  worth  noting  that 
shear  magnitudes  in  the  model  were  all  calculated  using  the  original  method  of 
deformation  averaging  (used  by  Kwok  et  al.  2008  as  discussed  in  Chapter  IV  Section  B), 
which  means  that  many  short-lived  shear  features  present  in  the  model  are  likely  filtered 
out  through  this  process.  Also,  with  the  overestimation  of  regional  velocity  vectors 
presented  in  Figures  12  through  15,  it  is  interesting  to  note  that  small-scale  defonnations 
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taken  from  six-day  averages  do  not  often  produce  shear  values  beyond  the  magnitude  of 
those  represented  by  the  RGPS  observations.  Further  analysis  will  be  required  to 
determine  the  significance  of  anomalous  large-scale  sea  ice  velocity  fields  and  if 
improving  the  velocity  estimations  alters  the  amount  of  small-scale  deformations 
produced  from  model  output. 

Phase  agreement  among  RGPS  data  and  POP-CICE  model  output  is  still 
somewhat  reduced  due  to  the  previously  discussed  problems  with  direct  comparison 
between  the  two  datasets.  As  mentioned  earlier,  the  POP-CICE  S3  model  run  still  uses 
re-analysis  of  atmospheric  forcing,  which  is  still  a  model  product  characterized  by 
smoothed  wind  fields  and  low  spatial  resolution  relative  to  the  sea  ice  model,  and  more 
importantly,  can  only  be  roughly  compared  to  the  actual  atmospheric  conditions  to  which 
sea  ice  response  is  recorded  in  RGPS  data.  Nonetheless,  reasonable  phase  correlation  can 
be  observed  in  the  relative  timing  of  various  shear  events,  particularly  in  regions  along 
coastal  boundaries  (regions  SI,  S2,  and  S5  in  Table  3).  Improvements  in  the  phase 
correlation  for  particular  regions  may  indicate  that  certain  combinations  of  sea  ice  and 
forcing  regimes  common  to  the  defined  regions  are  better  represented  in  the  model. 
Additional  analysis  is  required  to  determine  if  such  trends  exist,  but  such  investigation  is 
beyond  the  current  scope  of  this  work.  Remaining  discrepancies  in  the  magnitude  and 
phase  of  shear  deformations  between  POP-CICE  and  RGPS  will  continue  to  be  analyzed 
in  future  RACM  model  runs  once  model  output  overlapping  with  the  period  of  RGPS 
data  is  available.  The  presented  analyses  however,  provide  the  first  significant  step  in 
diagnosing  improvements  and  skill  in  simulating  regional  and  seasonal  deformations  by 
high-resolution  multi-category  sea  ice  models. 
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VI.  CONCLUSIONS  AND  FUTURE  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  analyses  of  results  presented  in  this  thesis  show  that  model  skill  in 
representing  sea  ice  kinematics  has  improved  as  a  direct  result  of  recent  upgrades  to  the 
POP-CICE  and  RACM  models  within  the  NPS  Arctic  Modeling  Effort.  This  assessment 
of  modeling  capabilities  is  by  no  means  comprehensive,  but  is  intended  to  evaluate  and 
where  possible  quantify  a  portion  of  the  model  improvements  related  to  the 
representation  of  sea  ice  thickness  distribution  that  has  been  made  in  CICE  model  and 
how  it  influences  the  representation  of  sea  ice  kinematics.  The  sea  ice  dynamical 
response  to  external  forcing  in  the  Arctic  is  one  of  the  most  critical  indicators  of  the 
Earth’s  response  to  the  changing  climate  that  has  been  observed  and  reported  in  recent 
decades.  With  ice  extent  continuously  on  the  decline,  and  ice  thickness  diminishing 
whether  through  export,  melting,  or  redistribution,  a  new  regime  of  sea  ice  and  the 
potential  for  seasonal  sea  ice  coverage  is  likely  on  the  horizon,  if  not  already  there. 
Therefore,  the  ability  of  numerical  models  to  represent  and  predict  sea  ice  response  to 
anticipated  changes  in  forcing  is  paramount  to  those  who  live  and  operate  in  the  region, 
as  well  as  to  those  who  need  to  model  and  predict  responses  to  global  climate  change  in 
general. 

Regional  high-resolution  models,  such  as  the  ones  analyzed  in  this  research,  are 
currently  not  practical  for  calculations  on  a  global  scale  due  to  the  computational 
constraints  of  current  computing  capabilities.  However,  the  ability  of  regional  models  to 
represent  small-scale  processes  and  detennine  their  relative  influence  on  larger-scale 
climate  is  important  for  advancements  in  earth  system  modeling.  The  efficiency  and 
predictive  capability  of  global  models  relies,  at  least  in  part,  on  understanding  and 
representation  of  processes  at  smaller  scales.  Sea  ice  displacement  and  subsequent 
deformation  patterns  depend  on  a  number  of  complex  factors,  which  need  to  be  isolated 
to  determine  their  relative  importance  and  contribution  to  ice  thickness  predictions.  In 
this  study,  we  have  attempted  to  move  this  effort  forward  through  analyses  of  sea  ice 

dynamics  and  the  ability  of  coupled  models  to  reproduce  observational  estimates  of 
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deformation  fields.  A  continuation  of  this  effort  with  output  from  subsequent  model 
runs,  in  concert  with  a  synthesis  of  model  thennodynamic  performance,  will  better 
position  the  scientific  community  to  make  determinations  that  may  help  improve  sea  ice 
prediction  capabilities  in  the  future. 

Results  from  the  NPS  coupled  models  of  POP-CICE  and  RACM  reveal  several 
key  indications  that  new  model  parameterizations  are  contributing  to  the  ability  of  models 
to  represent  kinematic  features  in  Arctic  sea  ice.  Multiple  category  sea  ice  thickness 
models  begin  to  address  the  issue  of  resolving  the  relatively  small-scale  dynamic  features 
that  develop  from  sea  ice  response  to  forcing  mechanisms,  which  influence  overall  ice 
thickness  distribution.  Increased  ice  motion  from  these  models,  while  overestimated  in 
several  regions,  show  that  new  model  parameterizations  allow  the  modeled  sea  ice  to 
have  an  appropriate  physical  response  to  anticipated  forcing,  but  the  estimated  magnitude 
needs  to  be  tuned.  Fully  coupled  regional  model  runs  (RACM)  show  improvements  to 
displacement  fields  in  this  initial  analysis  by  decreasing  the  dynamic  time  step,  which 
allows  for  more  plastic  behavior  of  the  sea  ice.  This  behavior,  along  with  additional 
feedbacks  from  sea  ice  strength  and  coupling  of  predicted  atmospheric  and  oceanic 
influence  may  lead  to  more  realistic  representation  of  kinematic  features,  which  can  then 
be  compared  to  the  thennodynamic  response  of  sea  ice  to  detennine  their  relative 
contribution  to  ice  thickness  distribution. 

Improvements  to  the  small-scale  deformations  represented  by  the  POP-CICE 
model  is  an  additional  indication  that  sea  ice  thickness  is  an  important  parameter  to 
conectly  model.  Sub-grid  representation  of  multiple  thickness  categories  in  this  model 
allow  defonnations  to  occur  based  on  the  weakest  category  present  in  a  grid  cell,  which 
has  led  to  increased  resolution  of  deformations  derived  by  the  model.  The  relatively  high 
frequency  of  model  output  (compared  to  RGPS  data)  has  revealed  additional  short-term 
deformations  present  in  model  output,  which  could  lead  to  improved  insight  into  sea  ice 
behavior.  Improved  magnitude  and  phase  agreement  of  shear  deformations  found  in  the 
POP-CICE  model  over  seasonal  analysis  show  promise  in  the  ability  of  such  models  to 
improve  the  representation  of  sea  ice  thickness  distribution  and  the  response  and  survival 
of  sea  ice  in  this  environment.  These  enhanced  defonnation  fields  diagnosed  to  greatly 
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improve  upon  the  comparison  to  RGPS  observations,  will  help  modelers  understand  the 
sensitivity  of  ice  parameterizations  that  are  required  to  physically  represent  in  future 
model  advancements  under  different  spatial  and  temporal  resolution  configurations. 

The  improved  representation  of  sea  ice  defonnations  found  in  this  study  show  that 
ice  thickness  and  forcing  parameterizations  are  integral  to  the  ability  of  regional  and 
global  models  to  represent  sea  ice  kinematics.  Deformation  patterns  observed  in  the 
Arctic  are  a  direct  response  to  the  dynamic  forcing  the  ice  cover  experiences,  and  as 
physical  changes  to  sea  ice  properties  are  imposed  under  changing  climate,  the  ability  of 
models  to  represent  these  features  will  lead  to  better  understanding  of  the  underlying 
processes.  The  methodology  displayed  here  will  serve  as  an  analysis  tool  to  be  used  on 
future  modeling  efforts  as  means  of  diagnosing  modifications  to  physical  and  temporal 
parameterizations,  and  help  quantify  the  importance  of  incorporating  kinematic  feature 
resolving  models  into  global  climate  modeling  efforts. 

How  significant  a  role  the  dynamical  response  of  sea  ice  plays  in  determining  the 
sea  ice  production  and  redistribution  of  ice  thickness  is  still  being  investigated,  at  the 
same  time  as  thermodynamic  changes  in  the  region  likely  have  a  large  impact  as  well. 
By  honing  in  on  model  capabilities  to  represent  sea  ice  dynamics,  future  studies  will  be 
able  to  better  attribute  the  changes  observed  in  the  ice  pack,  and  determine  the  relative 
contribution  of  dynamic  and  thermodynamic  patterns.  Atmospheric  and  oceanic 
feedback  processes  that  impact  sea  ice  are  challenging  to  quantify  and  attribute  direct 
cause  and  effect  parameters  for  sea  ice  modeling,  but  this  effort  must  continue,  especially 
at  high-resolution,  in  the  face  of  increased  activity  and  significant  climate  change 
observed  in  the  Arctic. 

The  political  environment  has  changed  in  the  recent  past,  and  a  focus  on  the 

interdependence  of  national  interests  on  a  global  economy  would  be  a  logical  framework 

from  which  future  operations  in  the  Arctic  should  proceed.  Increased  demand  for  natural 

resources  through  the  proliferation  of  industrialized  nations  poses  a  national  security  and 

environmental  challenge  through  increased  operations  in  the  Arctic,  but  should  also 

provide  an  opportunity  for  increased  research  due  to  higher  demand  for  expertise  in  the 

region,  as  well  as  the  potential  for  improved  international  cooperation  through  the  goal  of 
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finding  solutions  toward  a  common  interest  of  managing  our  changing  environment. 
Future  responses  to  climate  change  will  be  dependent  upon  our  ability  to  improve 
predictability  of  Arctic  climate  change  and  properly  prepare  for  the  challenges  that  lie 
ahead. 

B.  FUTURE  RECOMMENDATIONS 

The  ability  of  numerical  models  and  RGPS  analysis  to  accurately  represent  the 
dynamic  environmental  conditions  found  in  the  Arctic  is  paramount  to  the  goal  of 
predicting  future  responses  to  climate  change.  Observational  trends  indicate  that  sea  ice 
in  the  Arctic  is  entering  a  new  regime  leading  to  the  potential  for  seasonal  ice  cover  in  an 
ocean  formerly  obscured  by  multiyear  ice.  A  combination  of  model  simulation  and 
observational  validation  of  these  trends  will  be  required  to  determine  how  the  sea  ice  is 
responding  to  this  new  polar  environment.  Observational  techniques  have  been 
developed  recently  for  in-situ  instrumentation  and  data  processing  capabilities,  which 
will  help  to  verify  current  representation  and  prediction  capabilities  of  the  Arctic  models 
analyzed  in  this  study.  To  build  upon  current  capabilities  of  modeling  techniques 
to  represent  the  physical  characteristics  of  sea  ice  in  the  Arctic,  the  following 
recommendations  are  made  for  future  research  considerations. 

1.  Model  Data  Limitations 

This  study  was  limited  based  upon  availability  of  various  model  output  results, 
which  were  limited  to  runs  through  2002  (PIPS),  2004  (POP-CICE),  and  the  early  1990s 
(RACM).  Future  studies  could  use  this  analysis  as  a  template  for  evaluating  extended 
runs  in  terms  of  model  sea  ice  drift,  deformations,  and  ice  thickness  distribution,  as  well 
as  sensitivity  of  those  fields  to  sea  ice  parameterizations  and  different  forcing 
mechanisms.  Honing  the  skills  of  choosing  appropriate  model  time  steps  and  coupling 
parameterizations  that  allow  for  improved  realism  of  deformation  would  not  only  enable 
better  model  validation  against  satellite  observations,  but  also  help  to  detennine  the 
relative  importance  of  representing  dynamic  and  thermodynamic  processes. 

Improved  model  resolution  could  also  play  an  important  role  in  detennining 
which  parameters  need  to  be  represented  in  global  models.  Increased  spatial  resolution 
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might  capture  dynamic  processes  that  are  currently  underrepresented  in  the  9  km  models 
used  in  this  study.  This  is  not  a  trivial  or  straightforward  process  however  as  reducing 
the  grid  spacing  brings  along  additional  challenges  that  alter  current  model  assumptions. 
At  a  resolution  on  the  order  of  a  few  kilometers,  sea  ice  cannot  be  treated  like  a  fluid,  as 
individual,  larger  ice  floes  become  resolved  so  ice  moves  toward  more  granular  dynamic 
behavior.  Anisotropic  modeling  techniques,  may  help  indirectly  address  some  of  the 
challenges  of  representing  dynamic  sea  ice  behavior  and  better  handle  the  linear 
kinematic  features  discussed  in  this  study,  so  such  an  approach  should  be  further 
investigated  and  tested  in  high  resolution  sea  ice  models. 

Finally,  a  combination  of  high  and  low  resolutions  model  results  should  be 
analyzed  in  comparison  to  determine  how  important  physical  representation  of 
deformation  processes  are  to  the  overall  behavior  and  thickness  distribution  of  sea  ice  in 
the  Arctic.  Current  efforts  towards  improving  ice  dynamic  models  must  work  in  parallel 
to  achieve  a  better  understanding  of  the  thennodynamic  behavior  of  sea  ice  in  order  to 
determine  the  relative  contribution  of  each  to  sea  ice  distribution.  Regional  models  are 
capable  of  resolving  dynamic  and  thermodynamic  behaviors  at  a  scale  not  currently 
possible  in  global  models.  Therefore  a  thorough  analysis  of  the  relative  contribution  of 
small-scale  processes  on  the  overall  behavior  of  sea  ice  in  global  models  will  help  these 
models  incorporate  the  most  pertinent  parameters  as  computational  capabilities  allow. 
Regional  model  validation  in  this  manner  will  also  help  to  develop  usage  as  an 
operational  tool  for  defense  downscaling  applications  as  predictive  capabilities  are 
improved  through  a  better  understanding  of  the  underlying  processes  that  drive  sea  ice 
behavior  in  the  future 

2.  Satellite  Observations 

In  this  study,  RGPS  analysis  was  limited  to  data  for  the  years  1996  through  2001, 
although  additional  data  from  this  satellite  observational  tool  exists  through  2007.  Once 
model  output  is  released  through  this  time  frame,  continuation  of  this  analysis  should  be 
conducted  to  extend  the  length  of  this  dataset  as  a  means  of  validation  of  the  conclusions 
made  in  this  research,  and  to  detennine  if  additional  trends  are  observed.  While  the 
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RADARSAT-1  mission  extended  well  beyond  its  initial  five-year  charter  and  has 
provided  data  for  over  15  years,  several  new  satellite  observational  tools  have  been 
launched  (e.g.,  RADARSAT-2,  CryoSat-2,  IceBridge)  or  are  in  the  planning  phase  of 
development  (e.g.,  IceSat-2),  which  will  bring  about  several  improvements  to  the 
observational  capabilities  satellites  can  provide. 

The  ICESat  program,  which  uses  laser  altimetry  to  determine  sea  ice  freeboard  as 
a  means  of  estimating  ice  thickness  and  volume  distribution,  stopped  collecting  data  in 
2009  due  to  the  failure  of  the  laser  and  inability  to  repair  it  in  flight.  With  the  new 
ICESat-2  satellite  not  scheduled  for  launch  until  2016,  the  IceBridge  program  was 
developed  as  a  means  of  continuing  the  data  collection  in  the  interim  by  flying  aircraft 
over  specific  regions  to  make  similar  measurements.  While  limited  in  range  compared  to 
the  ICESat  coverage,  this  technique  maintains  the  capability  of  observing  thickness 
behavior  until  the  launch  of  the  improved  ICESat-2  satellite.  An  additional  platform  for 
observation  of  sea  ice  is  CryoSat-2,  which  uses  radar  altimetry  similar  to  RADARSAT-1, 
at  improved  resolution  to  estimate  ice  thickness  and  potentially  ice  motion  characteristics 
of  sea  ice.  This  next  generation  of  satellite  observations  of  Arctic  sea  ice  will  likely  bring 
about  new  resolution  to  sea  ice  observational  techniques,  including  basin  scale  estimates 
of  sea  ice  thickness  distribution  and  volume.  These  capabilities  should  be  exploited  to 
improve  modeling  capacity  to  represent  sea  ice  behavior  in  a  similar  manner  as  the 
analysis  provided  here. 

3.  Additional  Model  Parameterizations 

Several  parameterization  efforts  require  additional  attention  in  order  to  improve 
sea  ice  modeling  techniques.  Consideration  should  be  made  to  improve  the 
parameterization  and  prediction  of  snow  cover  in  the  Arctic  as  it  has  a  large  impact  on 
sea  ice  behavior. 

Incorporation  of  new  ridging  and  melt  pond  parameterizations  that  have  been 
developed  in  the  CICE  model  should  be  tested  in  high  resolution  model  configurations 
and  results  analyzed  to  determine  the  improvements  made  by  characterizing  these 


86 


important  features  of  sea  ice.  Thickness  validation  efforts  must  continue  in  order  to 
better  understand  the  three  dimensional  response  of  sea  ice  to  environmental  forcing. 

Atmospheric  forcing  has  been  constrained  in  the  modeling  efforts  by  limited 
resolution,  which  likely  contributes  to  biases  in  the  sea  ice  response  in  certain  regions. 
Whether  through  smoothed  and  low  resolution  forcing  prescribed  in  POP-CICE  from 
ECMWF  or  from  other  global  re-analyses,  the  atmospheric  forcing  present  in  current 
models  lags  in  its  ability  to  reproduce  the  true  observed  behavior  in  the  Arctic.  Realistic 
atmospheric  forcing  through  improved  re-analysis  and  higher  resolution  representation, 
such  as  predicted  atmosphere  in  RACM,  will  help  improve  modeling  efforts  to  resolve 
smaller  scale  forcing  mechanisms  that  likely  drive  sea  ice  response  to  the  changing 
climate. 

All  of  these  efforts  in  combination  will  help  not  only  the  modeling  community 
better  represent  the  Arctic  environment  and  its  response  to  climatic  forcing,  but  will  help 
to  detennine  the  strengths  and  weaknesses  of  different  parameterizations  in  future 
modeling  efforts. 
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